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ABSTRACT 


This paper considers a joint replenishment inventory problem with a continuous- 
review (S, c, s) policy for the backorder case with Poisson demands and constant 
procurement lead times. 

Whenever item i's inventory level hits s; (reorder point) or lower it triggers an 
order so as to raise item i's level to S; (order up point). At the same time any other 
item j] With inventory level at-or-below its can-order point C; is included in the 
replenishment. 

A Poisson demand model with a queueing description of the svstem’s operation is 
analysed, and comparisions are conducted for joint versus individual orders in the case 
of multi-item problems, where joint replenishment of several items may reduce -setup 


costs. 
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I. INTRODUCTION 


A. DESCRIPTION OF THE PROBLEM 

The purpose in constructing a mathematical model of an inventory system 1s to 
use it as an aid in developing a suitable operating doctrine for the system. The 
criterion most frequently used for selecting the operating doctrine is that of profit 
maximization or cost minimization. In some cases the task of determining the optimal 
operating doctrine is so difficult that it is either impossible or uneconomical to 
determine the optimal doctrine, and instead, one optimizes with respect to some subset 
of operating doctrines. Occasionally, the mathematical model may be so complicated 
that it is extremely difficult to do anything analytically. In such situations, simulations 
are used to study various operating doctrines. | 

We are concerned with finding optimal ordering policies for a multi-item 
inventory system in which the inventory position is under continuous review. [n this 
thesis, the times at which demands occur are assumed to be generated by independent 
Poisson processes with intensities h: (annual demand rate for item 1). We assume that 
all shortages are backordered and will be satisfied from the next shipment. 

We assume that a holding cost of H; per unit time is charged for each unit of 
item 1 in inventory. When an order is placed, a setup cost for an independent 
replenishment of item i is assumed to be C; = A + a;, where A is a fixed cost and a; is 
a cost which depends on the item ordered. The cost of a joint replenishment of both 


item iandjis A + a, + a. For instance, if item i triggers a replenishment, the setup 


cost of A + a; is applied ¥ item i, and a; to the jointly replenished item j . 

This type of cost structure is particularly appropriate when a group of items is 
ordered from the same supplier, or if a group of items uses the same means of 
transportation, and transportation costs do not increase proportionally when the 
quantity transported is increased. Thus, when a replenishment is made, there is a 
major fixed cost independent of the number of items involved and a variable cost that 
depends upon the number of items [Ref. 1: pp. 278-83]. 

Because of this cost structure, it 1s clear that there may be an opportunity to 


reduce total variable costs if one takes advantage of joint replenishment. 


B. SCOPE 

Our interest in this problem stems from the “random joint order policy” proposed 
by Balintfy [Ref. 2], and from an (S.c,s) policy treated by Silver [Ref. 3]. For the 
backorder case with Poisson demands and constant lead times, and the specified 
customer service levels, we will compare the individual order policy versus joint order 
policy. 

However, for selecting the control parameters of (S.c,s) in Joint Order Policy, we 
have considered slightly different procedures suggested by Schacck and Silver [Ref. 4], 
and Silver [Ref? 3]. 

Two programs were written in order to conduct the simulation experiments. The 
first program determines the control parameters of each item in the group for both 
independent and joint replenishment policies. A second program evaluates the model. 
Detailed user instructions and the program itself are contained in Appendix A and B, 
while examples of their use are presented in Chapter II. 

The individual order policy is demonstrated onlv for the purpose of comparison. 
For the sake of simplicity, we compare the joint order policy to the individual order 


policy where the classical lot size formula is used for each item. 


Il. GENERAL FORMULATION OF THE MODEL 


A. ASSUMPTIONS AND NOTATIONS 
I. Assumptions 
We are dealing with a group of items where the cost of replenishing two or 
more items at the same time 1s less than the total cost of replenishing the same number 
in separate individual replenishments. The assumptions of the model are: 


1. There is fixed cost A associated with each replenishment, and variable cost a, 
associated with each item involved in the replenishment. 


2. Demands for item i are Poisson distributed with parameter A.. 
There are no quantity discounts. 


4. Inventory holding costs are proportional to the average dollar value of 
inventory. | 


Ca 


Service level is defined in. one of two wavs: probabilitv. of no shortage per 
ep emSe gy cycle. and fraction of demand to be satisfied directly from the 
2. Basic Notation 
We use the following notation: 
A: fixed cost per replenishment. (independent of the number: of items and 
units involved) 
a; : fixed cost for including item i in a replenishment. 


1 
EC;: expected relevant annual inventory cost for item 1. 


C. : ordering cost pemandess((C) = Aysaeas) 


Py : desired customer service level, measured in terms of the probabilitv of no 
shortage per replenishment cvcle. 

P, : desired customer service level, measured in terms of the fraction of 
demand satisfied directly from stock on hand. 

L_ : replenishment lead time in years. 

r :imventory holding cost rate. 

v; : standard unit price of item 1. 

s: : reorder or must-order point of item 1. 

c: : can-order point of item 1. 

S: : order-up-to point of item 1. 


Q; : order quantity of item 1. 


Poet 


, | Poisson demand rate for item i, in pieces per year. 
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I, : average on hand inventory level of item i. 
n : number cf items in the replenishment group. 

Nj : expected number of replenishments which involve item i per year. 
NT; : expected number of replenishments triggered by item 1 per year. 


Other notation will be introduced as they needed. 


B. INVENTORY THEORY BACKGROUND 

Before considering the coordinated control problem it might be helpful to briefly 
discuss the related single item models whose solution will be a key element in 
determing the parameters of the (S, c, s) control system. 

1. Economic Order Quantity Model 

The E.O.Q. model 1s based on the assumption that the entire lot is added to 

stock at one time, and that the stock will be withdrawn at a constant rate, and no 
stockouts are permitted. 


Total annual variable cost = Ordering cost + Holding cost 


The annual order cost is obtained by multiplving the number of orders per 
year (A;/Q;) by the cost of placing an order Sak 

The average annual holding cost for item 1 is the average inventory (Q;/2) 
times the annual unit holding cost (H;= v;r). 

To obtain the minimum cost lot size (E.O.Q.), take the first derivative of total 
annual cost with respect to Q: and set it equal to zero. Solving this for Q;, we get the 
E.O.Q. formula: 


Once the economic order quantity is known, the expected number of orders 


placed during the year, n. can be determined: 


nj= Aj / Qi 
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2. Exact Model for the Backorders case with Poisson Demands and Constant Lead 
imes 


{f demands are random, the probability of stockcuts exists, and we must 
include a stockout ccst in the expression for expected annual variable costs. This 


expression is given in [Ref. 5: eqn. (4-61)] to be 


tenet Ela) B(O-sy 
Where 
@; : the cost of backorder per unit quantity of backorder. 
t, : the cost of backorder which is proportional to the length of backorder time. 
M, : expected lead time demand for item iin pieces, M;= A.L. 
H; : holding cost of item 1. 
E(Q;, s;) : expected number of backorders incurred per vear for item i. 
B(Q;, s;) : the steady state expected number of backorders for item 1. 

It is not easy to derive the optimal Q; and s; by use of the above exact 
expression for EC;. Fortunately, it turns out in practice that it is seldom necessary to 
use the exact formulation, except for the case when it costs very little to incur 
backorders. The approximation discussed in the next section has been shown to give 
good results when backorders are infrequent. 

3. Approximate Form with Known Backorder Cost 

The average annual variable cost equation can be approximated as follows 

[Ref. 6: pp. 59 - 62]: 


EC, = 4,C. Qo (Q:/ 2s: Vinson sy (2.4) 


where 


@; : the cost of backorder per unit quantity of backorder. 


E(M;> s;) :the expected number of backorders at the end of cycle. 


To minimize EC;,, take the partial derivative of the EC; with respect to Q; and 


s, and set these equal to Zero. 


a 4 ae) ey “ 
dEC;/dQ: = = AC, OF i Ah2 = (DA; E(M: > SQ; = 0) (2.5) 


i 


ee us: = i ra (@:A; / Qi)de(M: > os 0 


-—— 
L 
aN 

ee” 


The first equation(2.5) yields the optimal Q; for a given reorder point s:: 
wey 22(C + g E(M.> s.) / Hi] (2.7) 
The second equation(2.6) yields : 
dE(M:>s;)/ds. = - )) P. (M;) (2.3) 


where 
E(M,>s:) : the expected number of backorders at the end of cycle. 


- ) P, (M;) : complementary cumulative distribution. 


This is the complementary cumulative distribution of the Poisson random 
variable, X, evaluated at - Therefore, solving the second equation for S; in terms of 


Q; gives the maximum allowable probability of a shortage during a lead time : 
F’°(s:) = (H;Q:)i(p:A;) (2.9) 


To find the optimal pair (Q;, s;) that minimize EC;,, iterative procedures can be 
applied, such as described in [Ref. 6: p. 61]. Hower, we do not apply this kind of 
procedure to find out the values of control parameters in (S, c, s) policy, because the 
order size in the (S, c, s) policy is not fixed. 

In (S, c, s) policy, the order size is determined as an order-up-to point minus a 
reorder point for the item which triggers the replenishment action. Order sizes for 
jointly replenishable items are order-up-to points minus on-hand levels which are less 
than or equal to can-buy points but more than reorder points. 


4. Some Useful Formulas for Obtaining the ance Backorder Cost / Reorder 
Point /Safety Quantity in Fixed Order Size System 


In any practical situation, it 1s very difficult to determine accurately backorder 
costs. They can include such factors as loss of customers’ goodwill{i.e., in the future, 
he may take his business elsewhere), or in the military supply system, the cost of 
having some first line weapon svstem inoperative because of lack of parts. Other 
contributions to shortage cost can be somewhat easier to measure. However, these 


are usually small part of the total backorder costs [Ref. 5: p.18]. 
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We shail assume here that there is a backorder cost associated with each unit 
backordered. 

A compiete cost nunimization analysis of the optimal values of the model 
parameters can be made only if a cost can be assigned to each shortage. As an 
alternative to specifying the shortage cost, one can examine a given inventory policy 
and obtain the imputed cost of shortage within it [Ref. 7: p. 338]. Rather than 
explicitly costing backorders, we will specify a service level constraint and determine 
the imputed stockout cost. 

Whenever an organization uses a service level to establish a reorder point, it 
effectually establishs a stockout cost. Associated with a given service level is an 
imputed or implicit stockout cost. It is a simple matter to determine the imputed 
stockout cost for a given service level from the probability of a stockout [Ref. 8: p. 
162]. 

a. Probability of no Shortages per Replenishment Cycle 

A service level based on the frequency of service per order interval or 
replenishment cycle will indicate the probability of not running out of stock during the 
replenishment. This approach does not concern itself with how large the shortage 1s, 
but only with how often it can occur during the lead time (replenishment cycle). It is 
defined as the fraction of the replenishment cycles without depletion of stock: 

number of order periods with a stockout 
Service level fraction per cycle = 1 - ————————————————————————___—_ 


total number of order periods 
] - ARN 2 S:) 


el v= s:) = P (s:) 


= | - (service level fraction per cycle) 


The term P,(s;) is the stockout level fraction per order cycle or the 
probability of at least one stockout while awaiting a supplier's delivery. It is also a 
measure of the fraction of lead time periods during which the demand will exceed the 
reorder point. The magnitude of the stockout is ignored with this approach. 


The following formula is developed to compute the imputed backorder cost: 
PAs) Sat) (219) 


where 
P(s;): probability of stockout for item 1. 
H; : holding cost of item i, (H;= vjr). 
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b; : imputed backorder cost cf item 1. 


To obtain the appropriate stockout cost, solve eqn.(2.10) for b;. An 

example will illustrate the procedure in a later section. 
b. Fraction of Demands Satisfied Directly from the Shelf 

The fraction of units demanded and immediatelv filled from the shelf can be 
defined as the ratio: (number of units supplied)/(total number of units demanded) 
which is equal to {I - [(number of units short)/(total number of units demanded)]}. 

The above relationships must be measured over some fixed time period 
such as the order interval. 

For example, to obtain the stockout level fraction for units demanded: 
E(M; > s:} 


quantity demanded during a cvcle 


Stockout level fraction for units = 


Pees.) / Q. = (H.:Q:) / (b:A;) e)) 


where 
E(M: > s-): expected number of backorders for item i during the cycle. 
H; : holding cost of item 1, (H;= vir). 
b. : imputed backorder cost for item i. 
c. Safety Stock 
With backorders, there is no loss of sales since the customer waits for the 
arrival of the next order if stock is not available. The expected safety stock is defined 


as: 


Safety stock = }' (s - M) P,. (M) ee 
= s) P.(M)- ¥ M P,(M) = s- E(M) 
where 


E(M) : the expected lead time demand. 


O,if M <-s. 
ideas atyi. = S. 


The number of backorders per lead time 


i 


So, the safety stock 1s simplv the reorder point minus the average lead time 
demand. The reorder point is determined so that the customer service objective is 
satisfied. 


C. MATHEMATICAL STATEMENT OF THE DECISION PROBLEM FOR 
THE JOINT REPLENISHMENT MODEL 


1. Cost Equation 


The basic cost equation in the joint replenishment model is as follow [Ref. 3]: 


The first term represents the inventory holding costs and the other two 
represent the set up charges allocated to item 1. The total expected costs per vear of the 


group of items are given by 
= Vv — te a oa 
EC = Y EC.= ¥ (Lvir + NTA + Niaj) (2.14) 


2. Service Level Constraint 
The service constraint takes one of two forms depending upon the measure of service 
used. 
(1) Probability of no shortage per replenishment cycle 
This is simply the probability that demand for item 1 during the lead time will 


exceed the reorder pomm- 


P.(s) = P,(M;>s) = ¥ P.(M:) (2.15) 


Where 
P_(s): probability of stock out. 


Wl : lead time demand for item i. 


(2) Fraction of demand satisfied directly from the shelf 


The expected stockout quantity during the lead time 1s given bv 
E(M > s:) = )) (M; - s;) P. (M;) (2.16) 
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wnere 
E(M,>s;) : expected number cf stock out for item 1 in pieces 
during the lead time 


M, : jead time demand for item i 


A backorder occurs if and only if the lead time demand (which is Poisson with 


parameter M;) is greater than the inventory level of the item when the order is placed. 


D. INDIVIDUAL ORDER POLICY 

Under the the individual order policy, there is a fixed order quantity for item i 
that is ordered every time the reorder point is reached . 

Safety stock is needed to protect against stockouts between the time the reorder 
point is reached and the order received. As presented in [Ref. 3], the individual order 
policy (called independent control) is a special case of the joint replenishment policy (S, 
c, Ss) where Ge bye In this case item 1 is ordered only when the inventory position hits 


its reorder point, and the order quantity is always of size Q, = S.-s:. No other items 


+ 
will be included in the replenishment, hence a set up cost of A+a,; is incurred with 
each replenishment. 

As mentioned in the introduction, the main focus of this thesis is to compare 
joint replenishment with independent control to determine the cost savings that can be 
achieved. 

Let us first consider independent control. Then, c; = s., and the optimization 


problem is: 
Minimize oS (S; - Si ae 1)/ Be (S; = pally Niese (21) 


eA a)h; / (S; - S;) 


subject to 
P.(X Ss | AL) 2 Py (2.18) 
Or 
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where, X 15 the lead trme demand with a Poisson distribution having mean AGL 


The left hand side of eqn.(2.18): P(X S s;| A;L) 1§ the probability of no shortage 
per replenishment cycle. | ; 

In the leit hand side of eqn( 2.19 thie terra E(X> s;)'Q: is the stockout fraction 
of units demanded. Therefore 1 - E(X>s;)/Q; is the fraction of demand satisfied 
directly from the shelf. 

The exact solution to this problem requires that the pair (S;, s,) which mininuze 
eqn.(2.17) subject to either eqn.(2.18) or eqn.(2.19) be found for each item i. However, 
it have been observed by manv researchers (see: for example Silver and Wilson [Ref. 9]) 
that a much simplier sequential approach does almost as well from the standpoint of 
minimizing costs. This sequential approximation approach is what we describe. 

1. Applied Solution Procedure 

Stepi ss let Oss aaa 


the service level constraint. Then given this value of S; aS: the constraint is used to 


Find the value of Q: which minimizes EC; ignoring 


find the lowest allowable value of s;. This first step leads to the usual EOQ expression 


Step2 > Determine the reorder poms. 

(1)For the service measure, probability of no shortage per replenishment cvcle, 
find the smallest value of's;, which satisfies the condition of eqn.(2.18): 

P(X 3 sj AL) 2 Py 

(2)For the service measure, fraction of demands satisfied directly from the 
shelf, find the smallest value of 's;, which satisfies the following inequality: 

A-L-s:-ALP.(X 3 s.-1|A;L) + ssP.(X 3 s;[ AL) S Q{1 - P5) 


Steps) eu Ses 
Step4 : Compute the average annual variable costs. 


eh Ore een I) 2a ice A:L)vir Se ACor a:)A: / (S; - s;) using the values 


s, and S; estimated from Steps 2 and 3. A computer program for solving for these 


parameters 1s provided in Appendix A. 


Note that we can now obtain imputed stockout costs. 


Step5 : Compute the imputed stockout cost 
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(1)For the service measure, prcvabilitv of no shortage per repienishment cvcle: 
Be eeeani.(2-110): P-(s) = (HQ) / (b:A. wand solvetfor b,, be= H;Q,/ (1- Py) A: 

(2)For the service measure, fraction of demands satisfied directly from the 
shelf, it is necessary to solve for E(M > s;) and find the associated P.{(s) before using 
eqn.(2.11). It can be done alternatively by finding the complementary cumulative 
distribution of the lead time demand, introduced by eqn. (2.8). Then eqn.({2.11) vields 
b= H.Q./[- > P. (M;)] /A;. 


2. Numerical Examples 


Repeated from Silver's example in his paper [Ref. 3]. 


Example 1. Service Level = Probability of no shortage per replenishment 


cycle. 


Consider the following example involving four items: 
A = $50 . | 
a= Semronra = 1; 2.3.4 

= 2s VT. 


r 
L = | month (1e. 1/12 vear ) 


Ne iee | 
INPUT DATA 


item A; V; Py 
(pieced/vr) (S/pidce) (desired service level) 


6.90 0.95 

1.20 La 

BU 0.95 | 
20 0.95 | 





Use of eqn. (2.15), (2.17), (2.18) gives the results shown in Table 2. 
Example 2. Service Level = Fraction of demand satisfied directly from the shelf. 


Consider the same data as example 1, but with a new the measure of service 
level. Suppose the specified fraction of demands to be satisfied directly is 0.99. The 


results shown in Table 3. 
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TABLE? 
INDEPENDENT ORDER POLICY: SERVICE LEVEL 95% (P}) 





a) eet pet 
IIIDNn 
L.C~)=— 
nna 
CJ mt HCD 


TOTAL: 441.203 


TAD eae 
INDEPENDENT ORDER POLICY: SERVICE LEVEL 99% (P-4) 





E. JOINT ORDER POLICY 
1. Introduction 

Let us now consider the possibility of joint replenishments. As discussed 
earlier, we can reduce ordering setup costs by combining the orders of several items 
whenever an order is placed. This may result in increased holding costs but the savings 
in ordering costs may be greater than the increase in holding costs resulting in a net 
decrease in total costs. Further, it is likely that the occurance of stockouts will be 
reduced since we will frequently reorders some items before the inventorv levels hit the 
normal reorder points for those items. 

We do not seek to determine the “optimal” joint replenishment policy. 
Rather, we restrict attention to a simple, intuitively appealing, policy and seek to find 
the optimal values of parameters for the selected policy. The policy we consider is the 
(S, c, s), “can buy, must buy” policy which requires an order to be placed for item 1 up 
to the level S; any time the inventory position for item 1 reaches or falls below s;. In 
addition, however, we order any other item 1 whose inventory position 1s below the can 


buy point c,, we always order up to the level S.. 
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2. General Concepts for Obtaining the Values of Control Parameters 

Suppose a single item 1, is in a group of coordinated items and assume initially 
that the replenishment lead time is zero. The assumption of unit sized Poisson 
demands together with a zero lead time implies that s;= 0. 

Now item i, from time to time, will be faced with the opportunity of a 
replenishment at the reduced setup cost of only “a;". Such an event is caused by 
another item triggering a replenishment. Occasionally item i will trigger an order when 
its inventory hits the zero level. Under such circumstances we would like to determine 
the (S, c) pair which minimizes the expected costs per unit time. 

Silver shows that the opportunities to replenish at the reduced cost occur 
probabilistically according to a Poisson process with a rate [ft per vear, where Ul is the 
expected number of orders triggered per vear by all other items in the group. 

Let NI; be the number of orders triggered per year by itemi. Then ar is given 
by a » NT, where the summation is taken over alli * j. 

_ For example. if there are four items in the group and we are considering item 
2, we would have 5 = NT} = NT; +NT 4. 

It is very difficult to determine the value of ; algebraically. We will consider 
some heuristic methods to approximate this. 

Biever, {Rei 5: ean. 13] determines that the relevant cost equation {the 


subscripts i have been suppressed) is given by: 
EC.={S-c + p(l-p°)/(1-p)}! (S-cyxS ++ 1) vr/2 (2.20) 


+ p[c-p{l-p°)/(1-p)] vr/(1-p) + ApS A + da} 


where 


p= Ay (KK + I) 


The algorithm suggested by Silver in [Ref. 3] 1s a simultaneous optimization 


procedure to minimuze eqn. (2.17) with respect to S;,c:,s. The procedure is iterative. 


iG 
Solution procedures to determine the control parameters of the (S, c, s) policy 


are explained in the next chapter. 
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Hil, DEVELOPMENT OF THE MODEL 


In this chapter we consider the general joint replenishment problem with positive 
leadtimes. We obtain values of the parameters Se and S; using an iterative 
procedure recommended in [Ref. 3]. Silver [Ref. 10] derives the following expression for 
the total expected costs per year for item 1 (we suppress the item subscript for 


simplification of notation): 
ey SS Ee | 9°) / (1 - p)}7) {(S eS - c + 1) 2 (3.1) 


+ p[c-p(1- p*)/(1 - p)lvr/(1-p) + ApS A + da} 
. where . 


p= A(X + pW) 


p is the probability that a particular occurance is either a demand or an 


opportunity to replenish at reduced cost. 


A. SOLUTION PROCEDURE 
We recommend as a solution procedure a somewhat modified version of the 


procedure suggested by Silver. 


1. Algorithm 


Step 1: Iniualization 
Use the values computed by the independent control to initialize. Let k = 1, 


where k = index of iteration. 


Cmin ae Q. 

Cmax ~ Q 

EC(k) = EC(determined from the independent control) 
N(k) = 24 / Q™ 


H(k) = y IN, where allie jangle 


p(k) = A/(A + B) 
c(k) = 0 


ste 
\s 


S(k) = Q 
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Steas: Letk =k + 1 
c{k) = c{kK-1) + (c ea i 2a Ee(h) = EC(k-1) 
o(k) = c(k+1) - (Cray = Cin) / 2, if EC(k) > EC(k-1) 
S(k) = c(k)- p(1-py(1-p) + J {20 (a + Ad‘)/ver 
ecu (ernierp iepodimeempe yc =p)4} 
EC(k) = {p (1 - p®)/(1 - p)}"! fp [e - p (1 - pSV( - pvr 
(l-p)+ ApS A + ia} 


max — 


Sieees); Lest for convergence of c and S. 
Round the c and S values to the nearest integer. 
If c(k) = c(k-1) and S(k) = S(k-1), go to Step 4, else go to Step 2. And, 


repeat until S and c converge. 
meme. Repeat Steps |, 2 and 3 fomi = 1, 2, ..., n. 


Step 5:LetN, =ApS/{S-c+ p(l-p (1 - ph}. 


9 values are from the previous iteration and S and c values are from step 4. 
Step 6: Repeat steps 1-5 until Ic(k) - c({k-L)] <0 and |S(k) - S(k-1)| <0. 


Step 7 : Now determine the must-buy points, Ss for i = l, ..., n from the service 
constraints. 
(a) For Service Measure 1: probability (P,) of no shortage per cycle, find the 
smallest value s such that 
(1/p)P (X Ss+clAL) - P(X=s+ 1A L) 
- pst 1S P(X=xold L) (1/p)Xo = Pips 


(b) For Seas Measure 2: fraction (P5) of demand to be satisfied directly 
from shelf, find the smallest value s such that 
pee sgl Pa(X = $= IhX L) + s P(X SslAuh)} 
+ (1-p)°" SY p “0 (A L-wy-AL Pw, - Id L) 
+ Wy PAX SwW)} 
morc +p (i -p (1 -p)} (1 - Po) 


Step 8 : Recompute the values S; and c:. 


Sem. > tS: 
Oa Se 
Sop 
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2. Flowchart 





Initialization 


=“ C min 70 











ican a 
i 





Set trial of C 






Cleese 3 
C(jy=ec(k-1j aa 
2 






Compute S(k) and EC(k) 







(Use eqn.e at etep 2) 






Find pair of S,C;When lead time=0O 





(S$, ,C,) sin. EO(S,,¢,) 


Find Coordinated pair of S,C 


(Compute Ny ge ) 





¥O5EC (k-1) CEC (k) CEC (k- 1) 






Figure 3.1 Determination of coordinated pair (S,c). 
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TV. SIMULATION STUDY 


In order to determine how costs of ordering, holding, shortage, etc. are affected 


by joint replenishment, a simulation model was developed. The purpose of the 


simulation study was to test the cost savings by comparing the costs for the two 


inventory policies, independent and joint replenishment. The simulation was written in 


FORTRAN/77. 
A. ASSUMPTIONS FOR THE SIMULATION STUDY 

I. Poisson demands (exponential interarrival times). 

2. Items demanded one at a time. 

3. Constant lead times (independent of the size of the order and the number of 
items in a replenishment). 

4, For the joint replenishment case, place an order whenever any item hits its 
must-buv_point s: and order up to S-. Also. include item] * i if the inventory 
position for item j is below the can-buy point c.. Order up to Si. 

>. For, the independent control, order up to level S: whenever the stock on hand 
for item 1 reaches s-. : 

6. For a given item, there is never more than a single order outstanding in the 
joint replenishment case. 

B. MODEL INPUTS 

1. Number of items. 

2. Demand rate A; for each item. 

3. Holding cost rate v; for each item. 

4. Lead time (in months) L. 

5. Type of service level constraint (Typel, or Type2) 

6. Service level required. 

7. Unit price UP.. 

8. Group ordering cost A. 

9. Individual ordering cost a:. 

10. Stockout cost rate per unit stockout (percentage of unit price). 

ll. Time weighted stockout cost (percentage of unit price). 

12. S;, GS: for each item. 

13. Model options (Independent control or Joint Replenishment, etc.) 
14. Limits of simulation period, NSIM. 
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C. MODEL OUTPUTS 
1. Achieved service level for each item. 
2. Number of orders (independent and joint). 
3. Total ordering cost for each item. 
4. Total holding cost for each item. 
5. Total unit years of backorders. 
6. Total number of backorders. 
Total unit vears of stock held. 
8. Total stockout cost (fixed and time weighted) for each item. 
9. Total annual costs for each item. 
10. Total annual variable costs for all items. 
11. Standard deviation of total annual variable costs. 


12. A complete audit trail of all orders placed (optional). 


D. STRUCTURE OF THE PROGRAM 
l. Algorithm 


Step 1 : Read input data. process options and initialize variables. 


Step 2: Determine the earliest event time and type of event. 

Let TDA, be the time of the next demand for item 1 and TOA; be the next 
time of the for item of an order for item 1. 

Set the master clock time, CT, to the minimum of {TDA,, TOA}. 

If the event tvpe is a demand for item i then generate another demand time 
for that item. 

If the event type is an order arrival then increament the on hand inventorv for 


all items included in the order. 


Step 3+ Determine if anvonder is'seqiinear 
If the onhand level is less than or equal to the reorder point for that item and 
there are no outstanding orders for that item, place an order. If the ordering policy 1s 


independent, go to step 5. 


Step 4: Test for joint replenishments. 
If the onhand level is less than or equal to the can order point for another 


item ] in the replenishment set, place joint replenishment for the items to the level S.. 
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Step 5 : Determine if the simulation is complete. 
lint = NSU stom, ese Zo to siep Z. 


The program collects detailed statistics at each event time for unit years of 
stock held, number of stockouts, unit years of stockouts, number of orders, number of 
units ordered, and number of demands for each item. 

At the conclusion of the simulation additional summary statistics are 


computed for achieved service levels and costs. 


2/ 
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Figure 4.1 Simulation Flowchart. 
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a. V. EVALUATION OF THE SNIODE 


A. PERFORMANCE ANALYSIS 
1. Basic Examples Tested 
Table 4 shows the data used for the evaluation of the joint replenishment 
model. Several combinations (96 cases) of service level, service measure, and lead times 


were considered. The data include the test set considered in Silver’s paper [Ref. 3]. 


TABLE 4 
ITEM CHARACTERISTICS 


group number ITEM DATA 
_ | number of items A: (piece/yr.) — v: (S/piece) 
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2. Model Results 
a. Convergence of Parameter Values and Costs 
Table 5 shows the results of each iteration from Silver’s model for joint 


feplenisniment, Observe that five iterations were required for convergence. This 


fee, 


example considers service measure type 2 with a specified level 0.99. The lead time is 


one month. 


TABLE 5 
PARAMETERS AND COSTS 


Item | Iter. ee 
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b. Cost Savings: Joint Replenishment VS. Independent Control 
Table 6 shows the cost savings achieved by the joint replenishment policy 
for both the modified Silver's method and our heuristic method. 
The comparisions are provided for two cases: 
1. Service measure l, Py = 0.95 
2. Service measure 2, P5 = 0.99 
Observe that both the heuristic method and Silver’s method achieve cost 
Savings greater than 13 % by using joint replenishment, and the heuristic method 
results in greater cost savings than does Silver's method. 
3. Simulation Results 
Figure 5.1 shows the time history of stock on hand for items 1 to 4 generated 
by one simulation. The upper graph represents individual control and the lower graph 


represents joint replenishment. 
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TABLE 6 
COMPARISION OF TWO METHODS 


1. modified Silver's methed : P, = 0.95, L = 1 month 
item 0 P Cc S EC( JP) EC(IP) [Save(%) 


a cee ai ee ee a a nee 





1 0. 9973 32 119 vem | 21.6339 262.024 | 4.48 
2 oecacm | 4m 40 100 23. 261m Suen | nee, 
3 Oheyaa | 9° 52 97 70.9459 g8.860 | 20.16 
Pa o.9835 * 13 Sa \eScu| G70 172ml S24. 9771) 20.95 
EOtalmeece| 26omenc Ou 4416203) 13.19 
Peeaeuristic method : ee o> > eee month 
item , 5 C S EC( JP) EC(IP) |Save(% 
i sect si) | 110 177 2 oe 232024) sane: 
2 mlTS ISS 34 93. | 21.6035 Se Sg. || Benes 
3 Deocsoul 9 | Ge 94 Ben iis uC Omron oy 
mamno.9774.|, 14. 59 7S. hy WSS cl | Weir: 
total cost | 351.8467 | 441.203 | 16.83 


1tem p S C S Bea) Be GLE ) Save( %) 
1 peoogs | 25 |) 112 Ge, Ala. eye 220.984 | 4.08 

mm o.9467 | -1~ 35 95 ponOe71 see | eee 

3 0.9744 | 5. 4g 93 67. 82593 85.740 | 20.89 
Mito.983s./ 74 Ta SC) Abe oo 22 "| a) 82. 22a ty 21.be 
total cost | 366.2788 {| 423.323 | 13.48 


4 heuristic method : P, = 0.99, L = 1 month 
Peel p = C S) EC( JP) Be( TP) Save, 


1 0.9942 BS 104 ae 7a 202.9998 Ze02 984 8.14 
2 O92 7 2 -1 30 89 2076235 34.383 39, 26 
5 C27650 5 a2 Eo S555 912 85.740 Zoe5 6 
Bw HiO.9774 8 i So ea) 2G awe 62.61 23 I ee. Zig Yi 25.25 

Geicaecost | ool,0467 1 423.323 | 16.88 
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The simuiation results for several combinations of leadtime, service measure 
and sefvice level are represented in Tables 7 , 8,9. 

When lead time is less than 1 month, the heuristic method vields more savings 
in total cost and also satisfies the required service level. But, as the lead time becomes 
larger, it sometimes fails to meet the required service level, although it does continue to 
generate greater cost savings than Silver's method with the cost parameters examined 
in this thesis. Both methods generate substantial savings over independent control. 
In short lead times, both methods tend to overprotect against stockouts. They 
generally produce service levels in excess of what is required. 

One explanation for the better performance of the joint replenishment models 
for the short lead time examples is the fact that a kev step in the joint replenishment 
algorithm for determining the S. c, s parameters assumes that lead time is Zero. 
Adjustments are then made to those values to account for the nonzero lead times but. 
apparently those adjustments are not adequate. — Additional work is required to 
examine this issue and to explore other methods of adjustment. 

We should keep in mind the fact that the model satisfies all requirements in 
the experiment does not mean that it is a correct model -- merely that it is a plausible 


cone which has not been found inadequate by the data or experiment. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS a | os 

For each of several test combinations, we determined the independent and 
coordinated control parameter values and the associated expected costs. We simulated 
several years of operations with the derived joint replenishment ordering parameters to 
evaluate the effectiveness of the models as compared to independent control models. 
We also use the simulation to compare a simple heuristic joint replenishment model. 
With the simulation we were able to accumulate costs (ordering, holding and stockout) 
and effectiveness as measured by two service levels. 

In summary, the following statements can be made: 


1. Coordinated control results in total costs significantly lower than under 


. — 


independent control. The average cost savings are about [8.7 %. 

2. Coordinated control provides better service than independent control. Ina few 
cases, the actual service for an individual item is slightly lower than required: 
when this occurs, the item involved is the one that triggers the majority of the 
replenishments in the group. 

Cost savings of coordinated control increases as the ratio a/A decreas es. 

4. Cost savings improve as the number of items in replenishment set increases: the 
Mietemitems there are in the group, the more beneficial coordinated 
replenishment becomes. 

5. For a given service measure and lead time, the cost savings increase as the 
desired service level decreases. 

B. RECOMMENDATIONS 

The procedure to determine the parameters for coordinated control requires only 
the data for independent control. We observe that the joint replenishment models 
perform best for shorter lead times. This is probably the result of the fact that the 
joint replenishment algorithms initially assume zero lead times and then adjust the 
derived parameters heuristically to account for nonzero lead times. This adjustment 
procedure may be inadequate. Additional work is needed to consider other adjustment 


methods. 
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APPENDIX A 
“PARA”:PROGRAM TO SELECT PARAMETERS 


wes VARIABLE DEFINITION ¥** 


DM : ANNUAL DEMAND 
MR : REORDER POINT 
$/S$1 : ORDER|-= UP POLun 
C/Cl : CAN ORDER POINT 
IPC : OPTIMAL PAIR OF C 
IPS : OPTIMAL PAIR OF S$ 
: : ECNOMIC ORDER QUANTITY 
CI : AVERAGE ANNUAL TOTAL VARIABLE cost ODER) 
EC/EC1 : AVERAGE ANNUAL TOTAL VARIABLE COST(JOINT ORDER) 
TLM : LEAD TIME (IN MONTH 
IER : LEAD TIME(IN YEARS 
TL : AVERAGE LEAD TIME(IN YEAR) 
DLT : LEAD TIME DEMAND 
NIT : NO. OF TOTAL ITEM IN THE GROUP 
TN : NO. OF TRIGGERING AN ORDER BY ITEM I IN THE GROUP 
R / RHO : OPPORTUNITY TO REPLENISH AT REDUCED COST 
F : HOLDING COST RATE TO PROCUREMENT COST 
UP : UNIT PRICE 
OCI : INDIVIDUAL ORDERING COST 
OCG : GROUP ORDERING COST 
SVL : SERVICE LEVEL 
SCU : IMPUTED UNIT STOCKOUT COST GIVEN SVL. 
?SO : MAX. ALLOWABLE ROBABILITY OF STOCK OUT 
STON : EXPECTED NUMBER OF STOCKOUT DURING A ORDER CYCLE 
DIMENSION ITM(100) ,DM(100) ,UP(100) ,H(100) ,OCI(100) ,ECI(100) ,Q(100) 
Sc ae ,TN(100),$(100),1S(100) ,R(100) ,U(100) ,EC(100) ,MR(100). 
&10(100 DLT (100) ScU(100) .16(100) iO 
&C1(50,100),31(50,100) ,1C1(50, 100) ,1S1(50,100), IPC(100) , I1PS(100), 
&MR1(100),TL(100) 
* DATA INITIALIZATION 
STN=0. 
TVC=0. 
* READ PROCESS OPTIONS 
CALL OPTION(ID,LT,KEY,LTT,SVL,MSRCH) 
READ(ID,110)NIT,OCG,F 
Thiet io: 
DO 11 I=1,NIT 
Gt teat ae ,DM(I) ,UP(I) ,OCI(I) 
i 
ate .EQ. TL (Ts =TLT 
IF(LTT.EO.2) THEN 
CALL LEAD(TLT,TLY) 
TL(I)=TLY 
END IF 
CONTINUE 


* PRINT SELECTED OPTIONS 
WRITE(7,100)KEY,SVL,LTT,ID,MSRCH 
* CALCULATE 2AO=0 

ea) ak iL, 

DLT(1)= DMT) *TL (I) 
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G+ ect (I) 
S ey pes) OCI) 
INT (O(1 )+.5) 


* COMPUTE REORDER POINT:IND. ORDER POLICY 


1 : PRCS.OF NO SHORTAGE PER REPLENISHMENT CYCLE 
Z : FRACTION OF DEMAND TO BE SATISFIED DIRECTLY FROM SHELF 


CALL POISON(KEY,DLT(I),Q(1I),SVL,MX,M,MAX,CCDF1) 
MRI (T) =M 

IF(KEY.EQ.1) THEN 

Recut )= =e aeqh), (h.- SVL) “DM(1)) 


eke .BO.2) TH 
SCU(i)=H(2)#9(2) / (CCDFLADH(T)) 


END 
ECI(I)=H(1I)*(( +1.)/2.+(MR1(I)-DLT(I +OC*(DM(I Or! 
ECI(1)=t SOLVIT TS DET(G) LCL) O(L) tOC4) MRLCL) Bee ee). 
TVC=TVC+ECI (I) 

mii | Y=DM(1)/O(1) 

STN=STN+TIN(I 

CONTINUE 


Motte CieoOe ) IVC 


Oc=0c 
Q(1}= 
is )= 


te 


_* DETERMINE THE PARAMETERS FOR JOINT ORDER POLICY 


(REORDER / CAN ORDER / ORDER-UP POINT) 


Bee2 It=1 Na 
Weae STN- NCE) 


BO: 
R(I)=DM(I)/ (DM(I)+U(I)) 
* SEARCH LOOP FOR DETERMINIG A PAIR (S,C) : ITEM BY ITEM 


IT=1 

SEC=0. 

SEC1=0. 

DO 4 [=1,NIT 
pooch; tL)=R(T) 
K= 


CMIN=0. 
(Qo, 
C(K a 
S(K 
EC(K)= SR¢t(r) 
TEST=10.% 
RensSh=1 
LOOP=0 
K=K+1 
ieee .EQ.1)THE 
C{K)= e(R-1)4 (CHAR: CMIN)/2 


UF (TCaSE. .EQ.2)TH 
K)=C(K-1)- (CHA x- CMIN) /2 


END IF 
S(K) = C(K)-R(I)*(1-R(I)**C(K 1-R(I) )+SQRT(2*DM(I)*(OCI(I 

K) a nace IAACCRDS/ACE) HERG (R oR TES ARS, 1°R(r) 

-R(I)* pe) (TSG oe CARDTL)}/ (I-R I) )**2 
EC(K)=(1./(S(K)-C(K)+R(I)*(1.-R(1)* c(K))/ 1.-R(I ) 
* S(R)“C(R) )7(S(R) +C(R) $1 5H z}/eeR(T) - K)-R(I)*(1.-R(I) **C(R) ) 
/(1.-R(I)))*H(1)/(1.-R(1) )+DM(1) *R(I) **C(K) XOCG+DM (I) *OCI(T)) 
EeCRcana(e(S 
= +. 
IF(C(K).EQ.C Rei}? AND.S(K).EQ.S(K-1))GO TO 
Pe mece) GT. 0 9995%EC (K 1)).AND.EC(K).LE. EC(K- ie GO 106 
IF(EC(K).LT.EC(K-1)) TH 
MECeCK) Gl .cCK- 1) ) THEN 


Be, 


aa 


QAA A 


ar 


Os 


ZS 


CHIN=C(X) 
ICASE=? 
cleo NE 
ELSE 
CMAX=C(K) 
ICASE=2 
GC TO 5 
END IF 
ELSE 
IF(C(K).LE.C(K-1) ) THEN 
CMIN=C(K) 
ICASE=1 
GO TO 5 
ELSE 
CMAX=C (K) 
ICASE=2 
GO TO 5 
END IF 
END IF 
IF (LOOP .EQ.0)THEN 
CMIN=0 
CMAX=C (2) 
ICASE=2 
LOOP=1 
GO TO 5 
END IF 


* FIND OUT THE PAIR (S,¢), MiN@iMmZzeES EC 





DO 7 J=1l 
1F(EC(3), LE.TEST)THEN 
TEST=EC(J) 
ECt (7, eet sec(g) 
cU(t, Tae 
IT, Il)=SlJ 
aad IF 
CONTINUE 
LCL (T, I }=INT sii i+. y 
IS1(L, II )=INT fg ea 
SEC1= =SECI+ECI Hae. 
CONTIN 


SHAVE= (tve- wae 


WRITE(7, $07 )21 

DO 20 I=1,NI 

WRITE (7, 58) r, RHO(I,II),1C1(I,12),1S1(I, 11) ,EC1(1,11) 
CONTINUE 

WRITE(7,509)SEC1,SHAVE 


IF(MSRCH.EQ.2) GO TO 42 
COMPUTATIONS NEEDED ONLY FOR S. METHOD 1 
IF(II.GE. 2) THEN 
DO 21 I=1, 
IF (C1(T, tt). "EO. C1(I,II-1).AND.S1(1,11).FO.S1(1 11-1) con uamer 


IF(II.GE.2)THEN 
DO) 22 1) ee 


TROhheGr. NIT.OR. ECIGr ie a or 99957 he senor 
& 


p 
,EC1(I, 11) -LEJECI( Ebel mcommomaz 


END IF 
IF (I1.GE.7) GO To 42 
00 22 0t=) Nn 
TN(I) = Pie: Pe ; ~R(E)*4C1 (1 


(1) **CiCr re 
CONTINUE 
S 


Tipe Si(l 11) =2Ci( teh s RC 
)/ (1S RCE) $) 


40 


24 
ZS 
26 


42 


an 


MAAN 


28 


100 


0 
1 abal 
200 


502 
503 


504 
505 


DO 24 I=1.NIT 
SUN=SUM+TN(T) 
DO 25 Bel sMiT 
U(I)=SUM-TN(T) 
DC 25 T=1 NIT 
R(T)=DM(I$/(DM(I)+U(I)) 
LI=II+i 
SECI=0. 
GO TO 41 
WRITE(7,503) 
DO 27 I=1,NIT 
5 aoe 
St )=S1(I,Ii!l 
Tes (T}=1HT(s(1)+.8) 
IPS(I)=INT(S(I)+.5 
120g oe §04)1I,1PCc(I PES 1} 


WRITE(?. 200) 
* COMPUTE REORDER POINT:JOINT ORDER POLICY 


DO 28 I=1,NIT 
R(T) =RHO(I, IT) 
IF(KEY.EO. 1) THEN 
CALL HUSTI (R(T), ony, DLT (I), SVL,MBUY) 





MR(I)=MBU 
Se c(2) 4H (1) 
Soot LHR I 
PECs j=INT<C ae 2} 
Mo(ly=INT(S(i)+s5 
EG(E}=ECL Lt ea uit one DELEL AHL) 
END IF 
IF(KEY.EQ.2)THE 
GALL mUST2 (R(T), Ojo l),DLT(1) ,SVL,MBUY ) 
MR(I)=MBUY 
SS ee P 
5(1)=S(1I)+MR(I 
IC( I )=INT See 5 
IS(I)=INT(S(1I)+.5 
— PEE I)=EC1(1I,II)+(MR(I)-DLT(I) )*H(I) 
WRITE(7, Seat ROGER EL ewet fb), 1S(£),EC(L) 
SEC=SEC+EC(I) 
SHAVE=(TVC-SEC)/TVC 


WRITE(7,506)SEC, SHAVE 
SLOP 


FORMAT (1X, / Ps POBLEGMEDVOPLLONS 3 '7/,1X,' (* ee SVL. 


& =',F4.2, one Un eee a TYPE =',I2,2X,'* DATA='! ; 

&2X, ix Ss, METHD=' | Oe ae ee x -REY= 1 : SVL. BY PROB. OF NO SHOR 
&TAGE PER REPLENTSHMENT CYCLE! ee 2eeesovi. BY PRAGHEBOR DET 
© TO BE SATISFIED DIRECTLY FROM SHELF wie Gl T= 1. > CONSIANE 
&LEADTIME 2 : RANDOM (UNTFORM W/ MEAN ) ' /,' * DATA=1 : DATA 
&1 2 : DATA2 3 : DATA3 & DATA4 '|,/,' SEARCH METHOD =1 : 


& MODIFIED SILVER METHOD Be HEURISTIC! , 
Sid 20X,'*** INDEPENDENT ORDER POLICY RRR 

&,//,1%, ‘ITEM', Soule DEMAieimee 2, TBOO(N=INT.)',2X; “CINTEGER)', 
&2X, 'REORDER' , 2X, ‘EXP. BE an; '$/0 COSig) 

FORMAT J an, oe 2h) foc) 

HORMAG( ESnen,f7.3,3%,F5.2,2%,F 5.2) 

FORMAT (//, 20X, \kkk JOINT ORDER POLICY cE Rae alae 





SOK, 'REORDER', OG CANE UY © | aE Se eager eeeecleD Cosi 
sEORHAT(1X, ES: 4X, F7.3, ae FS. ny okies. 64,14,7%,253,4X, 59.3,3%, 


3) 

FORMAT (31%, 'SUM OF AVG. ANNUAL Gah F9.3) 
FORMAT('' //,15X,'** OPTIMAL PAIR sj ae 
S16X,'ITEM'.7X,'OPT. C'.7X,'OPT. ae 

Sel fs 13,9X,1I4,9X,14 

FORMAT (2X,13,4X,F6.4,7X,14,9X,14,7X,14,7X,F10.4) 


4] 


JIN 


FORMAT (45%, 'TOTAL COST s="—Fl0.4) Soman) 

FORMAT (! | ES OSG Ss PARAMETERS & EX.COSTS OF vOiINT REPLENTSHMENE! 
ae Pes (WHEN LEADTIME = 0)',/,15X,'* AT [GEPSDIONg Loess 

Six, 'ITEM', 4%,'eHO!',4X,'OPT. eu, De ‘OPT. S' 3K, 'EAPncesa.)) 

FORA ar one F6.4, eG ra 4X, 13 PN 4) 

SORMAT 23X, 'TOTAL =! LD Gag? A tS 

END 

SUBROUTINE POISON (KEY,DLT,Q,SVL,MX,M,MAX,CCDF1) 


PSO=1-SVL 

IF(KEY.EQ.1) THEN 

CALL SERV1(DLT,PSO,MAX,M) 

RETURN 
END IF 

IF(KEY.EQ.2) T 

CALL SERV2(0, DLT, PSO ,MAX,M,CCDF1) 
RETURN 


SD euS 
END 


SUBROUTINE MAXD(DLT,MAX) 
Dee a PMF (750) 


T= 
PE (1) EAP (- DLT ) 
Bo ae a ) 


GO 
BME (I) = =(DLE/(r~ 1))*PMF(I-1) 
CDF=CDF+P i) 

IF(CDF. GE. a Ege) GO TO 3 
I=I+l 

GO. ons 

MAX = I+l 

RETURN 

END 


SUBROUTINE SERV1(DLT,PSO,MAX,M) 


DIMENSION PUF (200), CDF (200) , CCDF(200) 
PMF (1)=EXP(-DLT) 
CDF(1)=PMF 1) 
CCDF(1)=1-CDF(1) 
CALL MAXD(DLT, MAX) 
DO 1 I=2,MAX+ 
ea re (DLT/(1- 1) )*PMF(I-1) 
CDF(I)=CDF(I-1)+PMF(Z) 
CCDF (I)=1-CDF(T) 
DO 2 I=1,MAX+1 
: IF (PSO. .LT.CCDF(I) .AND.PSO.GT.CCDF(I+1))GO TO 3 
RETURN 
END 


SUBROUTINE SERV2(Q,DLT,PSO,MAX,M,CCDF1) 
DIMENSION PMF (200) ,CDF (200), CCDF(200) 


STON=0*PSO 
ae =EXP(-DLT) 
CDE =PMF (1 


) 
CCDF(1)=1.-CDF(1) 
CALL MAXD(DLT, MAX) 
DO 1 I=2,MAX+1 
ie i: (DLT/(I-1))*PMF(I-1) 
CDF(I)=CDF(I-1)+PMF(I) 
CCDF(I)=1-CDF(I) 
CONTINUE 


T=0 
IF(I.EQ.0)THEN 


42 


aD 


~ sees 

FSr=DL Pavewereor( «)+1*CDF(I+1 ) 
END it 
pee - bas cON) ene N 


M= 
CcbF1= CCDF (M+1) 
ELSE 


I=I+l 
GO TO 2 
END IF 
RETURN 
END 


SUBROUTINE LEAD(TLT,TLY) 


CALL LRND(719325,U,1,1,0) 
Pivas2*TLT*U 

RETURN 

END 


SUBROUTINE MUST1(R,C,DLT,SVL,MBUY) 


DIMENSION PMF(400) ,CDF(400) 
IC=INT are: 

RC=(1./R)**IC 

RHS= SVL/(R**IC) 


BEAL )=EXP(-DLT) 

IF (PMF(1). LT.1.0E-70) PMF(1)=0. 
CDF(1)=PM F(1) 
DO 1 I=2,400 

PME(I)=(DLIT/ (I-1))*PMF(I-1) 

TF(PME(1L).LT.1.0E-70) PMF(I)=0. 

CDF(I)=CDF(I- 1)+PMF(1) 
CONTINUE 

JS=-2 
5 UM=0 
KF=JS+2 
KAL=JS+IC 


D0 2 IX=KF,KL 
SUM= RORY SRT eb xX) 
TEST=RC*CDF(KL4t1 “PHF (JS+2)- (R**(JS+1))*SUM 
IF(TEST.GE.RHS) THEN 

MBUY=JS 


RETURN 
END 


SUBROUTINE MUST2(R,C,S,DLT,SVL,MBUY) 


DIMENSION PMF(400) ,CDF(400) 
IC=INT(C+.5 
RHS=(S-C+R*(1-R**C)/(1-R))*(1-SVL) 
PMF(1)=EXP(-DLT 
IF(PMF(1).LT.1.0E-70) PMF(1)=0. 
EDEL )= BuE(1) 
Dol a—2 40 
pur(r) (DE / (2 1))*PMF(I-1 
F(PMF(1I).LT.1.0E- 70) PMF ( 
eae ee 1)+PMF (I) 
CONTINUE 
JIS=-1 
SUM=0. 
DO 2 IW=JS+1,JS+IC 
IF(IW.LT.1) THEN 


1)=0. 
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50 


Syl 
52 


35 


54 
5 
56 
37 


VcDS=0. 
ou iW.EQ.0)VCDF1L=CDF(1) 


VCDF=CDF (IW) 
VCDF1=CDF (IW+1) 
BN Dene 
SUM = Ut Greece ae - IW-DLT*VCDF+IW*VCDF1 ) 
TEeee ha ee ae 
V 


IF( JS.EQ.0)VCDF1=CDF (1) 
ELSE 


VCDF=CDF (JS) 
VCDF1=CDF (JS+1) 


PND to 
TEST=R**C* (DLT-JS-DLT*VCDF+JS*VCDF1 )+(1-R)*R**(C+IS ) *SUM 
IF(TEST. LE.RHS) THEN 

MB 


SUBROUTINE OPTION (ID,LT,KEY,LTT,SVL,MSRCH) 
Be OD 


READ(5,51) ID 
WRITE (6,53) 
READ(5,51) KEY 
WRITE (6 ,55) 
READ(5,57) SVL 
WRITE (6,52) 
READ(5,51) LT 
WRITE(6, 54) 
READ(5,51) LTT 
see 
READ(5,51) MSRCH 

RETURN 

FORMAT('1',' * ENTER OPTION FOR DATA FILE ( 1- 4) ! ',/,' _— 
&DATAlL 2 = DATA2 3 =DATA3 4 = DATA4') 

sh Gee | 

FORMAT('1!,' * ENTER THE LEADTIME IN MONTH ! ',/,5X, '1 = 1 MONTH 

',/,5%,'3 = 3 MONTH',/,5X,'6 = 6 MONTH') 

FORMAT('1',' * ENTER OPTION FOR THE METHOD OF SVL. MEASURE ! ',/,' 
& 1 = PROB. OF NO SHORTAGE PER REPLENISHMENT CYCLE ',/,4X,! 2 
S=_ FRACTION OF DEMAND TO BE SATISFIED DIRECTLY FROM SHELF ') 

FORMAT('1!,' * ENTER THE LEADTIME TYPE i) /,5X, '1 = CONsmam 


& ',/,5%,'2 = RANDOM (UNIFORM W/ MEAN) ') 

FORMAT('1',! * ENTER [He DeStese Se MEGS LEVEL ! a! 4k, ' (SHOULD 
&BE IN THE FORM OF "F4.2™" 2 Exp Ue Oooo LET 

FORMAT('1',' * ENTER OPTION | FOR. SEARCH METHOD (1- AYE vis 1 = 
&MODIFIED SILVER METHOD 2 = HEURISTIC') 

CO ee 


90.000 S20 Loses 
41.000 1.205 10200 
Ci arielele 2. SOs ORGS 
122.5000 PAPE EOrmoo 
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Sebnnce CFF 

&FN = PARA 

&FNO = &CONCAT OF &EFN OUTPUT 

&TYPE Do you need to compile your program ? (Y) 
SREAD VAR &R_COMPILE 

&IF &R_COMPILE NE Y &GOTO -RUN 

“-H FORTVS &FN 

&IF &RC EQ 0 &GOTO -RUN | 

Meee Our program did not compile; check for errors. 
&TYPE Do you wish to view the program LISTING file? (Y) 
&READ VAR &RSP1 

&IF &RSP1 EQ Y BROWSE &FN LISTING A 

&TYPE Do you wish to XEDIT the program file? (Y) 
&READ VAR &RESP1 

&IF &RESP1 NE Y &EXIT 1 

&COMMAND XEDIT &FN FORTRAN A | 

&TYPE Do you wish to run the program again? (Y) 
SREAD VAR &RESP2 

&IF &RESP2 EQ Y &GOTO -H 

Seatt 1 

-RUN 

FILEDEF O01 DISK &FN DATAI1 Al 

FILEDEF O02 DISK &FN DATAZ Al 

FILEDEF O03 DISK &FN DATA3 Al 

FILEDEF 04 DISK &FN DATA4 Al 

FILEDEF 07 DISK &FN OUTPUT Al (LRECL 133 

LOAD &FN (START 

&IF &RC EQ 0 SSKIP 9 

&TYPE Your program did not run Soy atc check for errors. 
&TYPE Do you wish to XEDIT the program file? (Y) 

SREAD VAR &RESP3 

poreeorkesrs NE Y &EXIT 2 

&COMMAND XEDIT &FN FORTRAN A 

&TYPE Do you wish to run the program again? (Y) 
&READ VAR &RESP4 

&IF &RESP4 EQ Y &GOTO -H 

S&EXIT 2 

eee LOUK OUTPUT IS IN THE FILE &FN OUTPUT A 
&TYPE Do you wish to BROWSE your output? (Y) 
&READ VAR &RESP 

&IF &RESP EQ Y &COMMAND BROWSE &FN OUTPUT A 
Steere Print your output file? (Y) 

&READ VAR &RESP7 

ae EQ Y &COMMAND PRINT &FN OUTPUT A 
&TYPE Do you wish to XEDIT the program file? (Y/N) 
&READ VAR &RESP5 

&IF &RESPS EQ Y XEDIT &FN FORTRAN A é 

&TYPE Do you wish to run the program again? (Y) 
&READ VAR &RESP6 

&RESP56 = &CONCAT OF &RESP5 &RESP6 

&IF &RESPS6 EQ YY &GOTO -H 

&IF &RESP6 EQ Y &GOTO -RUN 

onAlT 
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Boros need to compile your program ? (Y) 
$s FORTRAN COMPILER ENTERED. 11:58:02 
AAMAIN** END OF COMPILATION 1 ***x*xx 
AXPOISON** END OF COMPILATION 2 ***xxx 
ARKMAXD** END OF COMPILATION 3 ***x*xx 
AAXSERV1I** END OF COMPILATION 4 ***xxxx* 
AASERV2** END OF COMPILATION 5 ***%x%x 
ARLEAD** END OF COMPILATION 6 ***x*xxx 
AAMUST1** END OF COMPILATION 7 ***xx*x 
AAMUST2** END OF COMPILATION 8 ***xxx* 
AXOPTION** END OF COMPILATION 9 ***%*xx 
VS FORTRAN COMPILER CExr LS). sa exe ae 
EXECUITLTONSSEGiiiS 


* ENTER OPTION FOR DATA FILE ( 1l- 4) ! 
1 = DATA1 2 = DATA2 3 =DATA3 4 = DATA4 


* ENTER THE LEADTIME IN MONTH ! 
al 1 MONTH 


3 3 MONTH 
6 6 MONTH 


i 


* ENTER OPTION FOR THE METHOD OF SVL. MEASURE ! 
1 = PROB. OF NO SHORTAGE PERK REPLCENTSHE ie cre 


: 2 = FRACTION OF DEMAND TO BE SATISFIED SU URE CIE ene) shir in 
* ENTER THE LEADTIME TYPE : 
1 = CONSTANT 
2 = RANDOM (UNIFORM W/ MEAN) 


1 
* ENTER THE DESIRED SERVICE LEVEL ! 
soo ee BE IN THE FORM OF "F4.2" :; EX. ; 0.90 OR 0.99 ...ETC.) 
=) 
* ENTER OPTION FOR SEARCH METHOD (1-2)! 
1 = MODIFIED SILVER HETHOD Z = BEURT ore 


iL 
YOUR OUTPUP 2S iN THE FILE PARa Cuda 
Do you wish to BROWSE your output? (Y) 


Soll your output file? (Y) 

Do you wish to KXEDIT the program file? (Y/N) 
Do you wish to run the program again? (Y) 

R; T=1.47/2.39 11:58:39 
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xe SELECTED OPTIONS 


(* REY = 2 * SVL. =0.99 * LEAD TIME TYPE = 1 * DATA= 1 * S.METHD= 1 ) 


* KEY= 1 : SVL. BY PROB. OF NO SHORTAGE PER REPLENISHMENT CYCLE 
Zoe eo  aeeel  SORmDEMaAhDetOmbhwontvorthD DERECTLY FROM SHELF 

x LTT= 1 : CONSTANT LEADTIME 2 : RANDOM (UNIFORM W/ MEAN) 

* DATA=1 : DATAI1 2 3 Dawa 3 : DATA3 4 : DATA4 

* SEARCH METHOD =1 : MODIFIED SILVER METHOD Ze ee RES TIC 


xk INDEPENDENT ORDER POLICY KR 


eh LT DEMAND(LT. ECON.) (INTEGER) REORDER EXP. TVC($) ae CesT 
1S67506 hog Z2 220.984 a. 


Ih Zao O. eas eS 
2 Sol (Ome SS 143.178 143 3 34.383 11339 
a) 6.417(0.083 108.840 109 v 85.740 3499 
4 VOene 7 (0.083 eo oo 178 ae a4 Cal UP Ele 


SUM OF AVG. ANNUAL TVC($)= 423.323 


* PARAMETERS & EX.COSTS FOR JOINT REPLENISHMENT 
WHEN LEADTIME = Q) 
AT ITERATION : 1 
ITEM RHO ORIG OPT.S Bans Coon 
1 


0.9942 19 146 201.3498 
Zz ORIZ272 ea 90 Zi 7085 
& GEO 50 Zn) 85 66.0705 
+ 0.9774 45 bers: é 


e1e) 51218) S10 
TOTAL = 353.8584(0.1641) 


* PARAMETIRS & EX.COSTS FOR JOINT REPLENISHMENT 
(WHEN LEADTIME = 0) 
" ATAUTERATITON : 2 

ITEM RHO OPT.C OPT.S EXP.COST 


a G2397 2 84 he Om C7 
A 0.9499 36 a7 23724309 
3 0.9768 48 89 TO a 17 
4 Uagc50 a 145 he) o eiehebe: 


TOTAL = 370.8406(0.1240) 


* PARAMETERS & EX.COSTS FOR JOINT REPLENISHMENT 
(WHEN LEADTIME = 0 
x AT ITERATION : 3 
ITEM RHO OPT.C OPT.S  EXP.COST 


1 O99 7.3 84 eZ ZO iD oS 
2 0.9464 35 96 Zo, LOR 
3 0.9741 43 88 68.8164 
4 Oees5 70 143 Bo scors 


TOTAL = 368.3406(0.1299) 


* PARAMETERS & EX.COSTS FOR JOINT REPLENISHMENT 
WHEN LEADTIME = 0) 
AT ITERATION : 4 
in eyy ae OF tee ork. 5 Bar. COS L 


i 0.9974 87 eZ 210.8846 
Z 0.9469 36 25 23-1420 
5 0.9744 43 88 68.9584 
A Geee> 70 143 epote SlGle\s 


TOTAL = 368.8940(0.1286) 
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* FARAMETERS & EX.COSTS FOR JOINT REFLENISHMENT 
(WHEN LEADTIME = 0) 
7 he TL eRe ON ess) 


IT! RYO Co. © ore. S BAP COom 
1 Ceci 87 oe 2lOeeceD 
2 0.9467 36 96 coerce 
3 0.9744 43 88 68.9309 
t 0.9835 70 143 65.8689 
TOTAL = 368.7505(0.1289) 
xx OPTIMAL PAIR (C,S) ** 
ETEN SET. aC OE ar 
bt a0 Po 
2 BIS 96 
3 43 88 
4 19 143 
xkxX JOINT ORDER POLICY AKA 
TEEM RHO REORDER CANBUY ORDERUP EXPECTED ose 
1 OF 9ors Z2 ii2 Ve 21 Se 
Z 0.9467 Sul 30 25 22.0672 
c 0.9744 5 48 93 6758209 
4 0.9535 , : 1) Lo 


64.4122 
TOTAL COST = 366.2788(Omteaa, 
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[~ GenECTED OPTIONS : 
(* KEY = 1 * SVL. =0.95 * LEAD TIME TYPE = 1 * DATA= 1 * S.METHD= 1 ) 


* KEY= 1 : SVL. BY PROB. OF NO SHORTAGE PER REPLENISHMENT CYCLE 
2 eo ne er heel more DEMAND TO BE SATISFIED DIRECTLY FROM SHELF 
* LTT= 1 : CONSTANT LEADTIME 2 : RANDOM (UNIFORM W/ MEAN) 
* DATA=1 : DATA1 Zee DALAZ 3. 2 DATAS 4 : DATA4 
* SEARCH METHOD =1 : MODIFIED SILVER METHOD Ze Hib Un S.C 


*kX INDEPENDENT ORDER POLICY KKK 
LLIB LT DEMAND(LT. Se Se (INTEGER) REORDER EXP. TVC($) S/O COST 
1 24.167(0.083 156-000 sys. 33 232.024 1 


ec ib ebs' 
2 3.417(0.083 143.178 143 7 35.343 NOle 7 OZ 
3 6.417(0.083 108.840 109 ah 88.860 22.051 
+ FOete (0.083 34339 IL Rs: 16 84.977 13.453 


SUM OF AVG. ANNUAL TVC($)= 441.203 


* PARAMETERS & EX.COSTS FOR JOINT REPLENISHMENT 
(WHEN LEADTIME = Q) 

| * AT ITERATION : 1 

oa RHO OFLC OPT.S BAP soOon 


0.9942 is. 146 201.8498 
2 DegZiZ Si 30 a Aiea LOC = 
3 a O50 el 65 66-6905 
- 0.9774 45 iheys) O50 


030 
TOTAL = 353.38584(0.1980) 


* PARAMETERS & EX.COSTS FOR JOINT REPLENISHMENT 
WHEN LEADTIME = 0) 
BLS ITERAEION =:.2 
PET RHO OFT, Oris D205 = CO Rg 


Ih Oe 272 84 Poe 210.4767 
Z 0.9499 36 a7 23.4369 
3 0.9768 48 89 nOROs 7 
4 9650 at 145 66.88 


95 
TOTAL = 370.8406(0.1595) 


* PARAMETERS & EX.COSTS FOR JOINT REPLENISHMENT 
WHEN LEADTIME = 0) 
PIE tERATION =: 3 
TEM a Chee Orl.S 05 5 Bolas 


ih 9973 84 SZ ZI OR S550 
Z 0.9464 36 96 Zoe L011 
3 0.9741 43 88 68.8164 
* 0.9835 70 143 Go-oor 6 


TOTAL = 368.3406(0.1651) 


* PARAMETERS & EX.COSTS FOR JOINT REPLENISHMENT 
WHEN LEADTIME = 0) 
AT ITERATION : 4 
OEM RHO OFT GeOrT 5 ee Ose 


1 0.9974 87 2 210.8846 
Zz 0.9469 36 96 23.1429 
S 0.9744 43 88 68.9584 
4 We2o30 70 143 Bo. 9063 


TOTAL = 368.8940(0.1639) 
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* PARAMETERS & EX.COSTS FOR JOINT REPLENISHMENT 
{WHEN LEADTIME = 0) 
AT ITERATION : 5 


ITEM RHO OPT.C “ORime Ear COs. 
id 0297s 87 ley a0 feZs¢ 
Z 0.9467 36 96 29 wheat 
3 0.9744 43 88 68.9309 
4 039635 70 143 65n6 


689 
TOTAL = 368.7505(0.1642) 


xx OPTIMAL PAIR (C,S) ** 
M 


TEM OPT ae C20 ae 
1 87 ss 
Z S65 96 
3 43 88 
s 70 143 


xKX JOINT ORDER PORaGy Sate 


PTE RHO REORDER CANBUY ORDERUP EXPECTED €Osm 
i! 059973 a2 128 184 221 josse 
Z 0.9467 4 40 100 23.2078 
3 0.9744 5) 52 oy 70.9459 
4 059335 drs S3 sve 


67.1722 
TOTAL COST = 383.0190(03=aD 
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ANAANANANANAAANAANNANANANAANAAANANAANADAANANAANNANANANANANANAANANAANANANNYN 


APPENDIX B 
“SIM”: PROGRAM TO TEST THE MODEL 


VARIABLE DEFINITIONS 


honeee= OF LIONS FOR ORDER POLICY 

ieee OF LIONS FOR START INVENTORY LEVEL 
ferRoi— OPTIONS FOR PRINTOUT 

NIT = NO.OF ITEMS IN THE GROUP TO BE EVALUATED 
NSIti = SIMULATION, PERIOD 


Sey UG SIMULATION INTERVAL TIME (1YEAR) 

DM ANNUAL DEMAND RATE 

ME ORDER =UP EOUNT 

J) ECONOMIC ORDER QUANTITY (INTEGER) 

ne CAN ORDER POINT (CAN BUY) 

MR REORDER POINT(M1=IND. :M2=JOINT) 

TDA Pee oOreDEMAND OCCURING TIME 

TOA ORDER QOTY. ARRIVAL TIME 

M Mibe cs ORSEVENT FGR DEMAND OCCURED ITEM 

TA DIDEZ were ENT BORVOREDR OTYS°ARRIVED ITEM 

re ENDEA OF EVENT TYPE 

cr MASTER CLOCK TIME IN YRS(THE EARLIEST EVENT OCCURING TIME) 
DA Pete eenven SEIWEEN |Z SUCCESSIVE DEMANDS 

ARV UNIFORM RANDOM VARIABLE FOR LEADTIME WITH MEAN "ULT" 
ULT Boece leD LEADIIME 


i = ORDER QTY. 


Z TiDE PORMREPLENESH TRIGGER ITEM 

JT PIPE eenOne COUNT REPLENISHED ITEM 

IOH ONHAND INVENTORY 

rons San) LEVEL 

UOH UNIT YEARS ONHAND QTY 

UBO WEE sano SAcKORDER OTY 

NBO NO. OF BACKORDERS (QTY. ) 

NTBOQ TOTAL NUMBER OF BACKORDERS 

BOT BACKORDER TIME 

he Cr TOTAL BACKORDER TIME 

NOD ROA NUMeER OF ORDERS FOR EACH ITEM 
NDM TOTAL NUMBER OF DEMANDS FOR EACH ITEM 
NTDM TOTAL NUMBER OF DEMANDS DURING TMA& 
NJTO NUbEEr OF JOINT REPLENISHMENTS 

UP Ue eR LCE 

a HOLDING COST RATE 

HC HOLDING COST 

VC ORDERT& REVIEW COST 


OCI = INDIVIDUAL ORDERING COST 


eee JOINT (GROUP) ORDERING COST 

SCU BACKORDER COST PER BACKORDER UNIT 

SCT BACKORDER COST (TIME DEPENDENT FACTOR) 

wife AVG. ANNUAL TOTAL VARIABLE COSTS 

DIMENSION epee or 10) HR(1O) HC(10) HE (10), LOH(1Q).ASVL (19), 
&TDA(100) ,IT(10,100) ,JT(10,100),10(100) ,DA(100) , TOA(100 ,cEQ(10), 
2010) OCT (10), UF( 10} ee 0) .t 10) ae 1 OCd ee oO). 
&UOH(10,100),AU0H(100) ,UBO(10,100) ,AUBO(100) ,NOD(10,100) ,KB(10,100) 
&,ANOD(100) ,NDM(10,100) ,ANDM(100) ,NBQ(10,100) ,NTBO(10,100) ,KBT(10), 
&BOT(10,100) ,TBOT(10,100) , FRC(10,100),FRO(10,100) ,SVL(10,100), 

SAFR (100) 05 (10,100) -20C(100) COST (10,190) ,TCoST(100},NTOD Ooh 
S&HCOST(10,100),SCOST(10,100) ,SCTU(10,100),SCTT(10,100) ,NTDM(100), 
SAHCOST(100) ,ASCTU(100) , ASCTT(100) ,ABOT(10) ,ABOQ(100) ,ACOST(100), 


aAaAM 


TSGRT (100) OCT Plog Bilt 19) yea 
c(hO) , ‘TFRQ(100), NTDMS (100), NTODS 


COMO ICASE, 1SB@D 91, DAVaRV {ULL I, 1STP ,MRAME KEY 2EO 
ao = DATA INETIALT ASG ION 2oo-—— 


i. 
QO. 


0 
229375 


6) 
(100), 


TANDM 
TANOD 
TAOC 
Test 
STV 
ANJ 


----- READ PROCESS OPTIONS -----. 
CALL OPTION(ID,RATE1,RATE2,KEY,LTOPT,NSIM,ISTA, IP) 
----- READ INPUT DATA ---- | 


READ( ID ;HO2) NIT, OCG) EsveMsv aa 
ULT= eas 
DO 401 I=1, 
READ (ID, 100)DI(I) 
IF (KEY. EOF ane 
DE KEV or 
Scu 72 =RATE1*UP 
SCT(I )=RATE2*UP 
H{(I)=UP(1I)*F(TI) 
pea CONTINUE 


WRITE(7,198)ID,KEY,NSIM, ISTA, IP LIOR] (Ms eae en 


DOw i tai) 
1SEED(t)= x TSD+rAe3 


= 
tat 


OOOOOO000 0r O00 


2(I),MC(1), ME(1) , 1EQ(1) ,OC1(1) (UP( 1 rem 


~ 

to 

Hj 

~~ 
Hkhe 
——. 
OO0O00O 








10.**10 
on 


=~ 


Monn nn i 


Orj 
© 


DO a2 
DO 2 


I Ree HHH 


-~ 


Hc 
ao 
Ow 


HoOrrHoa 
CWANKO 
DOoOdOMN 


Za Mo 
Ke: GC 
Ow © 


| 
OOOO OODOODOO0000 
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( 
M 
M 
14 
( 
BOT( 
T 
Pet 
ae 
AFR 
( 
si 





CCl 


eee” 


Cc, C,> 
~~ we Cyne CG, 
~— 


pies HA DODDDVDDOVDWDODODVODODODODODDVODOOOCVOOOOCO00 000 
~—- TJ a ® e e ° e e e e ® e ® rt e 


DO sat 


© Buca —~q~— cq. 


i ote 


ores 


D.KEY.EQ.1)THEN 
Q(T) 


D.KEY.EQ.2) THEN 


tH 

| 

ti 
—~W 
tj 

He 

{1} 
WO 


ae 
TE 
ls 
M 
EZ 


~ -P- 


IF(ISTA.EQ.2) THEN 
CALL START 
IOH(I)=ISTP 

END IF 
3 CONTINUE 


goo PRINGOULTSINEUL DATA <-<-<-- 


D 4 I=1, 
4 ey 4) imolnen ee rhOCl kt?) Me Cl) ME(1),IOH(1),ULT 


AMA 


WRITE(7,97 
DO 5 I=1,NIT 
WRITE (7,98) 1,UP(I), F(T) ,H(Z) ,OCG,OCT (I) ,SCU(Z) ,SCT(I) 
= + 
TF _(IC.GE. 65.AND.IP.EQ.1.OR.IC.GE.65.AND.IP.EQ.2) THEN 


Meebo e (7,200) 


END IF 
CONTINUE 
= Sa e SLaAkheolULALLON ---- 
IF (IP.EQ.1.0OR.IP.EQ.2)THEN 
WRITE(7,2 3505 
: BND IF 


DO 777 ITER=1,NSIM 
IF (ITER.GT.1.AND.IP.EQ.1.0R.ITER.GT.1.AND.IP.EQ.2) THEN 
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noe _9 


END 2 


IMAXST MAX+S LMT 


e 
20 IF(CT.GT.TMAX) GO TO 77 
N = Ntl 
Z JOINT = 0 
«ee=- DETERMINE THE INITIAL DEMAND OCCURANCE TIMES FOR 2s Chews oss 


FE EE.) Re 
DORMS l=t ar 


ICASE=I 
CALL DEMAND 
ee =DA(ICASE) 
6 UOH(1I, ITER)=U0OH(1,ITER)+IOH(1)*DA(ICASE) 
c 
C ---- DETERMINE THE EARLIEST EVENT ----- 
C (DEMMAND OCCURED / ORDER PLACED / ORDERED QTY. ARRIVED) 
C 
1 (TON ah .TDA(2)) THEN 
ELSE 
TD=TDA(2) 
END IF 
DO 7 I=1,NIT 
TF (TDA(IS. LE.TD) THEN 
TD=TDA A(L) 
M=I 
END IF 
Ey CONTINUE 
LET=0 
CT=TDA(M) 
IOH(M)= TOH(M)=1 
IF(MSV.EO.1.AND.IOH(M). zQ. -1) KB(M, ITER)=KB(M, ITER)+1 
- NDM(M, iTERSNDM CM ITER) +1 
C ~--- SET NEXT DEMAND OCCURING TIME FOR THE ITEM :'M' ---- 


ICASE=M 
CALL DEMAND 
TDA(M)=CT+DA(ICASE) 
GO One 

END IF 


LPP LE.TDA(2)) THEN 
TD=TDA(1) 


TD=TDA(2) 
END IF 
DO 8 I=1,NIT 
1F(TDA(T) LE.TD) THEN 
TD=TDA (I) 


BNGeae 
8 CONTINUE 


IF(TOA(1).LE.TOA(2)) THEN 
TA=TOA(1) 


S 
. poo) 
DO 3" t=1, NI 
TR( Ron hy. "UE TA) THEN 
TA=TOA(I) 
IA=I 


END Ie 
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e 
c 


CONTINUE 
IF(TA.LE.TD) THEN 
GO TOMS 


sud 


Coa, 
ET=TD. 
ICASE=M 
CALL DEMAND 
cea ae 

IF(MSV.EQO.1.AND.IOH(M). Be -1) KB(M,ITER)=KB(M,ITER)+1 
Os ITER)=NDM(M, ITER) +1 


END 


M)=IOH(M)-1 


DETERMINE SHOULD ORDER ===-=- 


ee .GT.MR(M) ) THEN 


JOH 
REeue 20. 1) THEN 
WRITE (7, SOicl Lene uOn UM tLER)  (LOR( 1), [=1,4) 


M, ITER)=UOH(M, ITER) +IOH(M)*DA(M) 


END 
TECIP. EQ.2)THEN 
WRITE(7,601)CT,IET,M,UOH(M, ITER) , IOH(M) 
END IF 

ieee + 1 


1d (IC. Soo awe PEO.) .On LC .GE. 65. AND oie 506 2) een oN 


ie a0 
WRITE (7,200) 
END IF 


Goero 20 
END IF 


IF(K.EQ.1) THEN 
LET=1 


K=K+1 
IT(M,K)=M 
Teor) LT.0). THEN 


NBQ(M, ITER)= - IT0H(M) 


IF(KEY.EO. 4 10 i) = 03 Dk )+NBO(M, ITER) 
IF(KEY.EQ.2) IO(M) = ME(M)+NBO(M, ITER) 
Per toOrt. no. cay = CT + ULT 

IF( DIOP... £02 THEN 

CALL ARRIVE 

TOA(M)=CT+ARV 
END IF 


BOT(M,ITER)=TOA(M) - CT 
IF(NBO(M, ITER) .EQ.1)TBOT(M, ITER) = TROL Ms ITER)+BOT(M, ITER) 
eae ITER)= UBO(H, ITER)+BOT(M, ETER) 


IF(IP.EQ.1)THEN 


WRITE(7, gee) CT Pe, VERY Sako POOUCMe Ll reR)  UBOUMeerreR), 


IO(M) , TOA(M) , (IOH(I) ,I=1, 4) 


TF(IP. e972 .2)THEN 
WRITE(7 ee a ny M,NBQ(M, SER BOT(M,ITER),UBO(M,ITER), 
), TOA(M pels teat 
ENDSLE 
RE =—1C 


+ ] 
IF (IC. CeO oAND.LP-b@my.OR.LC.GE.65.AND.IP.EQO.2) THEN 


Ic = 0 
WRITE (7,200) 
END IF 


ELSE 


E(REY:£0.2) TOU) 


Th eTOPT . EO. ck 

IF(LTOPT.EQ.2)THEN 
CALL ARRIVE 
TOA(M)=CT+ARV 


= IE 
saa ee “1OH(i) 
Y=CT+UL 


5 


NAA Nn 


10 


END I 
fF (1P.EO.1) THEN 
i nRTEEC?, 503)CT,IET,M,IQ(M) , TORMM) , ( f@H{ 1) 2-2) 


iF(IB. fC .2) THEN 
ee 7,603)€T, TET ,M,1Q(M), TOAWM) , LOHGD 


"1 = 1c 
hei Ge ah "65. AND.IP.EQ.1.0R.1C.GE.65.AND.IPJEOg2 = Tae 


Le 
WRITE(7,200) 
ENDS 


PND ue 
NOD( MeITER) = "NODC 1 eR +1 
Ceii: ITER) = “OC(M, ITER) +OCG+OCI (M) 


BALLS 5 SHOl -3)60 Too 
IF(KEY.E0.2)GO TO 92 


ENDa aE 


~-=- DON’T ORDER UNTIL PREVIGUS” ORDER ARR eee 


TE(M. EO Tih) The! 
IET=0 


TFC IORI LTO) Bie 

See ITER)=" 10H (11) 

BOT(M,ITER)=TOA(M)-CT 

IF(NBO(M, ITER) .£0.1)TBOT(M, ITER )=TBOTC(Mp Tee eo iER 
UBO(M, pe =UBO(H, ITER)+BOT(M, ITER) 

IP (UR Os ly) PHET| 
WRITE(7, S04 er lia ily ay , ITER) , BOT (Me ITER), UBO(M, ITER), 


(IOH(I),I=1 
END IF 


IF(IP.EQ.2)THEN 
WRITE (7, 604)CT,IET,M,NBO(M, ITER), BOT(M, ITER) ,UBO(M, ITER) , IOH(M) 
END 
1C =IC+ 
IF (IC. GE. 65. AND.IP.£Q.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 


ne 
Re 200) 
END IF 


ELSE 
UOH(M, , TIER) =UOH (M, ,ITER)+IOH(M)*DA(M) 
IF(IP.EQ.1)THEN 
WRITE(7,501)CT, IET,M,UOH(M, ITER) , (IOH(I),I=1,4) 
END IF 
IF(IP.EQ.2) THEN 
WRITE (7, 601)CT,IET,M,UOH(M, ITER), IOH(M) 
END 


1 ie 
Me (I0.GE.68.AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 
ir = 


WRITE(7, 200) 
END IF 
END IF 


GO TO 20 
Bose 


DO 10 J=1,K 
sean ee neces) ofl < /Se TOSSzZ 
ET= 


K=K+1 
ee 

IF(IOH(M).LT.O)THEN 
NBO(M,ITER) = ~IOH(M) 
IF oe M) = IEO(M)+NBO(M, ITER) 
IF(KEY.EQ.2 = ME (it) +NBO (4, ITER) 
IF(LTOPT.EQ. 1 STOR \ = Coa 
IF(LTOPT.EO.2) THEN. 
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CALL ARRiVE 
an ee CT+ARV 
BOT(M ITER) = rey) =GE 


TE (NRO (M,ITER).EQ.1)TBOT(M, ITER)=T80T(M, ITER)+BOT(M, ITER) 
UBO (Me LTER)= =UB0(M, ITER) +BOT(M, ITER) 


Tene EO. .1)THE 
WRITE ( /502)CP, IET,M, NEO (M LTR), B /BOT (M ITER) ,UBO(M,ITER), 
3 — IO(M) ,TOA(M), (IOH(T 4) 


IF(IP.EQ.2)THEN 
WRITE (7 pe ay My Te a BOT(M, ITER) ,UBO(M,ITER), 
& ), TOA(M IOH(M) 
END IF 
re = reer 
ia (IC. a Soo ND eles hOnl On... Gh .Oo AND. IP.nOn2) THEN 


Ic = 
WAETE (7, 200) 


a 
OR CN oO). AR ie =IEQ(M) 
TPA 9, 64 3 5 OM =ME (M)-IOH(M) 


Mao TOrl. EO. 1) TOR (H)= =CT+ULT 
Pec rOrl.bO.2) TEN 

CALL ARRIVE 

mee CAAT) = CT+ARV 


LF(IP. P £0.) THEN 
wes) Chalbt-M1O(M) TOA(M) .( 10H(1)7I=1,4) 

END TF 
Bere ace ) THEN 

Prous GEATET,M,IO(M), TOA(M) , IOH(M) 
END IF 

Ze = IC 
)e (IC. oo i“ AND te. BOwrOkslC.GhiGoeAND.IP.EO.2)- THEN 


Lc 
WRITE(7, 200) 
D IF 
ENDS E 


NOD(M,ITER) = NOD(M,ITER)+1 
OC(M, ITER) = Ci, ITER)+OCG+OCI(M) 


END IF 
IF(KEY.EQ.1)GO TO 20 
: ---- DETERMINE THE JOINTLY REPLENISHABLE ITEM ---- 


91 DO il I=1,NIT . 
IF (TOH(T) GT.MR(Z) .AND.IOH(I).LE.MC(1) -AND.19Q(I).£Q.0) THEN 


IO(I) = SRB. 
IF(LTOPT EQ ew T) Sead optae 
IF(LTOPT. 

CALL ARRI E 

TOA(I)= CT+ARV 

END IF 
NOD(I,ITER)= NOD(I,ITER)+1 
OC(I,ITER) = OC(I,ITER)+0CI(I) 
JOINT = i 

iaalee TOR ace nee ) THEN 
WRITE (7, eee IT, Ton ac }) , ors aay 


Teeaerc 
IF ce GE. "68. AND.IP.EQ.1.0R.IC.GE.65.AND.IP.EQ.2) THEN 


IG. 0 
WRITE (7, 200) 
END 
END IF 
Hot CONTINUE 
NJTO¢(ITER) = NJTO(ITER) +JOINT 


Be 


GO TO 20 

92  IF(IOH(M).LT.O) THEN 

NBO(M, ITER)=-IOH(M) 

UBO(M, ITER )=UBO(M ITER) +BOT(M, ITER) 

BOT(M, ITER)=TOA(M 
TE (NBO(M, ITER) EQ. I STBOT(M, ITER)=TBOT(M, ITER)+BOT(M, ITER) 

WRITE(7,504)CT, IET,M NBO (M, ITER) ,BOT(M,ITER) , UBO(M,ITER), 

Rg (IOH(IS , I=1,4) 


END IF 
IF(IP.EQ.2)THEN 
WRITE (7, -604)CT, IET,M,NBQ(M, ITER) , BOT(M, ITER) ,UBO(M, ITER) , IOH(M) 
E 
new TCs 
IF (IC. cE. 65. AND.IP.EQ.1.OR.IC.GE.65.AND.IP.EQ.2) THEN 
WRITE(7, 200) 
ND IF 


ELSE 
UOH(M, TEER) suOH CH, ,ITER)+IOH(M)*DA(M) 
IF(IP.EQ.1)THEN 
WRITE (7, 501)CT, IET,M,UOH(M, ITER), (IOH(I) , I=1,4) 


END 

TF(IP, EQ.2)THEN 

WRITE (7, '601)CT,IET,M, wee ITER), a 
END 


IC We 1 
Ig (IC. ages “65. AND. IP.EQ.1.0R.IC.GE.65.ANDmnr B@e2 (STEN 
ue 
EE 200) 
NDEe 


C 
7 eee= SET NEW ON HAND, LEVEL POR THESE Se eee 


93 CT=TA 
IET=2 
TOH(TA)= IOH(IA)+IQ(IA) 
DO I=1,K 
ra. EQ. IT(TA, I) ) THEN 
IA.I)=0 
10 IA}=0 
NTBO (IA, ITER)=NTBQ(IA, ITER) +NBQ(IA, ITER) 
BOT(IA, ITER)=0. 
NBO(IA. ITER)=0 
TOR (TA)=10. *10 


ND 
LF (IA. EO. ee I) ) THEN 
JT s K)=0 


=0) 
Takin (IA, ITER)=NTBO(IA, ITER) +NBQ(IA, ITER) 
BOT(IA, ITER)=0. 

IA. ITER)=0 

IA)=10.**10 


MW CONTINUE 
ie sss ee ) THEN 
WRITE (7, ,500)CT ,IET, IA, (ICH@Ri-1 4 


E 

IF( Te .20. 2) THEN 

WRITE (7, 600)CT LED, laOn ae) 
END It 


J Le 1 
1G)9 (Ee a "65. AND.IP.EQ.1.OR.IC.GE.65: AND @gP O42 moHEN 
Lg 
WRITE (T, 200) 


END =e 
GO] TOmze 
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C 
2 ew--- COMPUTE INTER-SUMMARY (ANNUAL) ---- 


77 %F(IP.NE.4)THEN 
WRITE(7,300)ZTER 
END IF 
DO 5 line NI 
SE Teg 
NTOD(ITER)=NTOD(ITER)+NOD 
TCLTER “OCT (TIER) +OC(T, I 
TFS. FQ.1) 
oF (RB TER) tae 0) THE 
FRC(L,ITER)=KB(I, ITER) #1. 0/NOD(I, ITER) 
SE 


EL 
FRC(I,ITER)=0. 
N 
mer [tee = (1) -keeG@n ITER aa00. 


ations ,ITER).GT.0)THE 
_ERO( ITER )=NTBO(T, ITER) AL. O/NDM(I, ITER) 


“FRQ(T, , LTER)=0. 
END IF 


SVL(I, ITER)=(1.-FRQ(I,ITER))*100. 
END IF 
15 CONTINUE 
_ IF(IP.NE.4)THEN 
Hope i=1, NIT 
WRITE(7,507)I,NDM(I,ITER) ,NOD(I,ITER), vaeAE . ,ITER) ,UBO(1,ITER), 


seta 
GER 
INR, 


SNTBO(I, ITER), TBOT(I, ITER), SVL(I, ITER), IOH(I 
16 CONTINUE 

en oe Ea »NTOD( ITER) 

eeweee? 209 )NITOCITER 
phe 1S 


c 
GeeoeeecoieULE COSTS ---- 


DO 17 I=1,NIT 
HCOST(I, ITER) =UO0H(I, ITER)*H(I) 
scrt(t, TTER )=UBO (7 -ITER)* sch +H(I)) 
SCTU(I, ITER)=NTBO(I,ITER)*SCU(I 
SCOST(I, ITER) =SCTU(I, ,TTER) +SCTTAT TT ITER) 
COST(I, ITER) =HCOST (I. ITER)+SCOST(I, ITER)+0C(1L, ITER) 
17 TCOST (ITER )=TCOST (ITER )+COST(I, ITER) 
IF(IP.NE.4) THEN 
WRITE(7, S14) _ 
DO 18 I=1, 
18 WRITE(7, 615) I HOST (I ote) ScTU(1 
gy GC Mee AcOSMiT ITER) 
WRITE (7, Big voenleaene TCOST(ITER) 
N 
DO 19 I=1, 
19 KBT(I)= cRBT(I)+KB(T, ITER) 
77 CONTINUE 


C ---- COMPUTE FINAL-RESULTS (MEAN AND STD. DEVIATIONS) ---- 


WRITE(7,510)NSIM 
DOg@z. I=1,NIT 
BOr2 a Jat IER 1 


RaSreorisenr 1, 1lER), 


NEoDS(fy=NzOeS(T}<OE ET 
NTODS(I)=NTODS(I)+NOD(I,J 
TUOH(I)=TUOH(I SEE 
TUBO(I)=TUBO(I)+UBO(I,J 
sae: I )=TFRT(I peo ae. 

QT }=TERO I)+NTBO(I,J)*1.0 
oct I)=OCTT(I)+0C(T,J 
THCOST (I )=THCOST(1I)+HCOST(I, J) 
PoGiUGMy=TSeCTT( 1 )+SCIT(1I,J) 


BN, 


NN) 
a" 


Ze 
as 


Ze 


26 


Zi 


97 


98 

100 
8 
eke 


IDSs, 
200 








TSCTU(I)=TSCTU(I)+SCTU(I,J) 

[OC 22 252 NIT 

ANDE (2 )=NEDES (1) / 71a aX 

ANOD(TS=NTODS (1) /TMAX 

AUOH I(=TUOH I /THAX 

AUBO(I S=TUBO(T) /TMEX 

ABOT 1}= TER? I / THAX 

A800 IS=TE RO 1§/TMAX 

ANB(I) =KBT(1)*1.0/NTODS(I) 
AFRO(I)=TFR (1) /NTDMS I) 

TF (HSv.EO 1 asl (1 E ie + ANB(T))*100.. 
IF(MSV.EQ.2) ASVL(I)=(1.-AFRO(I))*100. 
AOC(L)=OCTT(1)/THAX 
AHCOST(I)=THCOST (I) /TMAX 

ASCTU(1 =ISCTU(T)/THAX 
ASCTT(I)=TSCTT(I {THA 

ACOST( 1 )=AHCOST(1)+ASCTU(I)+ASCTT(I)+A0C(I) 
TVC=TVC+ACOST(I) 

DO 23 T=1 ITER-1 

NTJITO=NTITO+NITO(I) 


IF(NTJTO.GT.O)ANJ = NTJTO/TMAX 
DO 24 I=] 1TER=1 
Stpveasont st) (23 EVG)2 2 
STDVC=SORT(STV/ (ITER-2) ) 
DO Zoe i=i. 
Ee r342)LaNDM(2), ,ANOD(I) ,AUOH(I) ,AUBO(I) ,ABOQ(I) ,ABOT(I), 
oo 26 I=1,NIT 
Ne 
TANOD=TANOD+ANOD (I 
TAOC=TAOC+A0C(TI) 
CONTINUE 





WRITE(7,512)TANDM, TANOD 
ae 513)ANg 
WRITE(7,517 
DO 627 Tait 
WRITE(7,515) I,AHCOST(I),ASCTU(I) ,ASCTT(I) ,AOC(I) ,ACOST(T) 
Were ace 
WRITE(7,518)STDVC 
STOP 


FORMAT(2X,' ', ///,' * INPUT DATA FOR COSTS:',/ 
& PRICE’ ,2X,'H.RATE!' ,2X,'H.COST',3X,' ORDER COST 
&/,36X, 'GROUP' ,6X, "IND. AR ALOTY tonya a area 

FORMAT (2X,12,4X,7(F6.2, 3X) 

FORMAT(F7 .3,2X,5(13, 2% R= 3(56.2" — 

FORMAT (I3,3X,F6.2,3X%,12, 3X, 2(11,3% ) 

FORMAT bf SK, ikk INPUT DATA :',/,3X,/,12K,'* USED DATA FILE NO. =! 
&,1I3,/,12X,'* ORDER POLICY OPTION =' I3,/,12X,'* SIMIMULATION LIMIT 
a YRS) =" 73. /,12%,'* STARE Leven OPTION =''13,/,12X,'* PRINTOUT 
SOPTION ='.1I3,/,12X%,'* LEAD TIME TYPE ="',13,//12%,///,7%) '=seRDER 
&POLICY OPTIONS' ,//,10X,'1 : INDIVIDUAL ORDER POLICY',/,10X,'2 : JO 
&INT ORDER POLICY',//,7X,'* START LEVEL OPTIONS' ,//,10X%,'1 : REGULA 
&R (IND. POLICY W/ E.0.0 ; JOINT POLICY W/ ORDER-UP QTY.)!',/,10X,' 
&2 : RANDOM (REORDER POINT + 1; ORDER-UP POINT)',//,7%,'* PRINTOUT 
SOPTIONS',//,10X,'1 : PRINT ALL ITEM',/,10X,'2 : PRINT ONLY EVENT 
SITEM',/,10X,'3 : PRINT ONLY SUMMARY'./.10X,'4 : PRINT ONLY FINAL R 
SESULT',//,7X%,'* LEAD TIME TYPES ',//,10X,'1 : COSTANT LEADTIME' ,/, 
&10X,'2 : RANDOM LEADTIME (UNIFORM W/ MEAN)',//, 6x, | Ux {SERVICE LEVEL 


hes eT SHO 2a, UNG 
x, 'B/O- COST, 


& > ie at ee claw ,oX,'* AVG. gEBeD TIME : YRS.) 5 ae 
Sl 2 MON.)!, ae 3X, ‘rte 3X, 'YR DEMAND', ate i ust But! | 
&3X,'CAN BUY! , a, 'ORD. UP TY! , Sk, oes CEVEL: Ae  EkP. Ine) 


ae CRY 12,04, 8/755 On, 4 ee SPO 3k, 13,1 3) 

FORIIAT ye 35X, 1kx SIMULATICN RESULTS aI 11,208, (EVENT TYPE 
& QO : DEMAND OCCURED ee ORDER PLACED 2: ORDER OTY ARRIVED) ' 
&,// Sha EVE ore ee 2x, 'BACK ORDER! ,3X,'UNIT YRS(QTY)' ,4x, 
&ON ORDER’ , 10K, (ON HAND! ,9X, ‘JOINT REPLENISHMENT ' i a TEMES 28, 


&' *TZPE', LOR 'OTY.',2X,'TIME' ,3X, 'ON HAND! ,2X,'B/O. S05 TY. i 
one (ax, ou ,4%,'2',4K%,'3' ,4K,"4" , 2h, ee 2, ‘ob: “or 8 
AR/T. 
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AGT 


‘eke 


Sl 
SZ 
Shs! 
514 
515 


S17 


518 


FORMAT('1!, fil 20%, pa INTER - Sat (ANNUAL ) ae ics 


&'aT THE YEA 6D) 3,//,3X, ‘ITEM #',3X,'# DEMANDS | , 1H Saee el 
© 10x," UNIT YEAR’ i rk, SBACKORDER | fa hea Ef of |, 40%, 
&'ON HAND',2¥,' B/O- 3X, "TIME! 5x, TEVEL(S) 


*CORMES | ' FG, 4,28, 12’ ae i r4aeX. 2 14, ry 


SORMAT(!' §'86.4,2X,12,3X,12,54X,14) 
FORMAT(' ' | F6.4,2%,12,3X,12,18X,F8.2,22Xx, 4(T4, T),) 
FORMAT(' ''F6.4,2X,12,3X,12,18X,F8.2,28X 
FORMAT(' ' F6.4,2X,12,3X,12,5X,13,3X, F6. BS x. F6.2,3X,13,1X, 
& F6.3,1X,4(14,1X 
oe Dee We 2. ok 125k, 1s 3k, fo.3,L11X,F6.2,3X%,13,1X, 
Agpier(’ § F6.4 2X,12.3%.12,37K,13,1%,F6.3,1%,4(14,1X)) 
Popter(' ' F6.4,2X,12,3X,12,37%,13,1K,F6.3,7K,14 
FORMAT ( ! ye 4°2X.12.3X,12,5K,13,3X,F6.3,11X,F6.2,14X,4(14,1X)) 
Bertier’ F6.3.2%,12,3%,12,5%,13,3%,F6.3,11X.F6.2,20X,14) 


PORraAl eo4,12,50%,13,9X,£6.3 


Owe eee 4 7X, 13,9X,F9.4,2X%,F6.3,;/2%,13,3X%,F6.3,5X,F6.2, 
& 


? 





FORMAT(1X,' TOTAL =' ,7X,15,5X,15) : 

FORMAT(1X,' TOTAL NUMBER OF JOINT REPLENISHMENT =! 75) 

FORMAT ('1! //1, 20X,'** FINAL apeULe (AVGS) 10K, 'DURING' ,I3,2X, 
S'YEARS  *** | 0 // 3X, ‘ITEM ts X,'# DEMANDS',3X,'# ORDERS' ,9X, ‘AVG 
&. UNIT YEAR! ,4X, AVG. al SERVICE ae nip 
&'ON HAND! ,3X,'B/O igh. 'TIME' 3X, 'LEVEL','(','%!,')*) 
FORMAT (SK. 12,6X, a x ey F ee Sal Bie, F7. eK home 2X, nae <o 

F6 
FORMAT(1X,'TOTAL AVG.=',F7.2,5X, F6.2) 
FORMAT(1X,'AVG. NUMBER OF JOINT REPLENISHMENT =' ,F6.3) 





BOniotlas,///,10K,' ~* COMPUTED ANNUAL COSTS (S)', 


Lee 
&' ITEM! ,3X,'HOLDING' ,9X, 'BACKORDER' ,9X,'ORDER' ,6X, 'TOTAL' ,/,22X,'QT 
oe 7k, ‘TIME! ) | 


FORMAT 
FORMAT 


Ze, SZ OP) «3 , 4X) ) 
B/S oe, WTOTAL = = reac s).G is, 





FORMAT 1X,///,10%,' * COMPUTED AVG. ANNUAL COSTS ($ ie e,@ 
by abe , 3X, eye Poe SE O@KORDER fore ORDER’, 6X, 'TOTAL? / ,22X,50T 
&Y TIME | 

FORMAT (1X, / Pook) SlenD. DEVLATION =" 2 Fs6ce) 

N 


SUBROUTINE DEMAND 


DIMENSION ISEED(10) ,DM(10) ,DA(100) ,MR(10) ,ME(10) , 1EQ(10) 
COMMON ICASE,ISEED,DM,DA, ARV, ULT,I,ISTP,MR,ME,RKEY, IEO 
CALL LRND ( TSEED( ICASE) U,1,1,0) 

_DA( ICASE)= -(1 ee SS one) 


&ND 


SUBROUTINE ARRIVE 


DIMENSION ISEED(10), HSA ARV OLE. TSSP MAH RY TES 
COMMON ICASE,ISEED,DM,DA,ARV,ULT,I,ISTP,MR,ME, KEY, IEQ 
CALL LRND(ISEED(T), hu iO) 

ARV= 2.*ULT*U 

RETURN 

END 


SUBROUTINE START 


DIMENSION ISEED(10),DM(10),DA(100),MR(10),ME(10), ee 
COMMON ICASE, ISEzD, DM,DA,ARV,ULT,1,ISTP,MR,ME,KEY 
IF(KEY.EQ.1) THEN 
ChE Hgts) CARE) 
IS= (IEO(1I)-(MR(I Gane 
ISTP= IS+MR(1)+1 
RETURN 


19) 
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30 
31 


Sz 


S32 
34 
35 
S65 


ay 
So 


40 


004 
230 
O41 
O77 
eZ 


END IF 
TE( Kez £9.20 


CALL Ean (2SEED(1) Lao) 
Is= (ME (T= (UR 1)+i$)*u 
ISTP= IS+MR(z)+1 

RETURN 

END IF 
END 


SUBROUTINE OPTION(ID,RATE1,RATE2,KEY,LTOPT,NSIM,ISTA, IP) 
RES io) 


Se eae 
WRITE(6,3 
READ (9, ieee 
WRITE(6,31 
READ (5, $4) KEY 
WRITE(6, 
READ(S, $4) LOFT 
WRITE(6, 
READ(S, 50) sim 
I ; 
ReaD(>,$4) Isms 
WRITE(6, 
ee ea IP 
RETURN 
FORMAT (I : 
Sea * ENTER OPTION FOR ORDER POLICY (1 - 2) t ' 7 Sie 
&1 = INDIVIDUAL ORDER POLICY' '/, 5X,'2 = JOINT ORDER POLICY '$ 
FORMAT('1',' * ENTER OPTION FOR START LEVEL (1 - 2) | 1) 7/05 


os REGULAR (IND.ORDER POLICY 2 ESCoOmy SUCTi] Bea eener Orne : ORDER U 


P OTY.)',/,5X,'2 = RANDOM BETWEEN (MUST ORDER LEVEL+1) AND (ORDER 
SUP LEVEL)" 


FORMAT('1',' * ENTER OPTION FOR PRINTOUT(1 - 4) ! ',/,5X, = PRI 
SNT ALL ITEM',/,SK,'2 = PRINT OUR MSE ali ee 1 8x, 3 = PRINT ONL 


&Y SUMMARY ',/,5X,'4 = PRINT ONLY FINAL RESULT! 3 


FORMAT (11) 
FORMAT( 1's = sete OPT ION GOR DATA PEs (1 = ie : A 

&5X,'1 = DATAL' ,\2 = DATAZ ses DATA3' ,3%, = path 
FORMAT('1",! e PYITER SIMULATION "LIMIT YEARS ri Ay eo ‘Ol = SIMU 


&LATION LIMIT IS 1 YEAR! ,/, 54, O02 Se cee re) LIMIT IS 2 YEARS! 


a n,m eS '99 = SIMULATION LIMIT IS 99 YEARS ! 

FORMAT ( | lie z’ ENTER THE “TYPE OF LEADPIME(1 =) 2) !',7,5%,'1 =sear) 
&STANT LEADTIME | ,/,5%,'2 = RANDOM (UNIFORM WITH MEANS" 

tk ee x ENTER RATE OF STOCKOUT COST TO UNIT PRICE Ul / pe 


Pein nom 8 dl = sTOchoun coon ee UNIT IS 125% OF UNIT PRICE’ ,/ aw 


& STOCKOUT COST PER TIME IS 92 [eOr SUIT Ear.) 
SO SZ abK)) 
h 


050.00 eh) 2 ik 
-000 025 025 112 177 %I59 O10600) 006. 30> UGrrZ 
UOC OOS -1 035 095 143 O10:500 O0l220 000s 
000 007 008 048 093 109 010.00 002.70) ecu 
.000 010 007 O77 I50 178 5010-00 Se07zs0 Socerze 
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&TRACE OFF 

&FN = SIM 

Seno — &CONCAT OF &FN OUTPUT 

&TYPE Do you need to compile your program ? (Y) 
&READ VAR &R_COMPILE 

&IF &R_COMPILE NE Y &GOTO -RUN 

=H FORTVS &FN 

&IF &RC EQ 0 &GOTO -RUN 

Saeeeee Our program did not Compile; check for errors. 
&TYPE Do you wish to view the program LISTING file? (Y) 
&READ VAR &RSPl 

&IF &RSPl EQ Y BSROWSE &FN LISTING A 

&TYPE Do you wish to KEDIT the program file? (Y) 
&READ VAR &RESP1 

peter e orl NE Y &EXIT 1 

&COMMAND XEDIT &FN FORTRAN A . 

&TYPE Do you wish to run the program again? (Y) 
&READ VAR &RESP2 

&IF SRESP2 EQ Y &GOTO -H 

&EXIT 1 

-RUN 

FILEDEF O01 DISK &FN DATAI Al 

PILEDEF O02 DISK &FN DATAZ Al 

FILEDEF 03 DISK &FN DATA3 Al 

FILEDEF 04 DISK &FN DATA4 Al 

PeEBDEF O7 DISK &FN OUTPUT Al (LRECL 133 

LOAD &FN (START 

&IF &RC EQ 0 &SKIP 9 

Biers Your program did not run correctly; check for errors. 
&TYPE Do you wish to XEDIT the program rile? (Y) 
&READ VAR &RESP3 

PbeeoReEoPS NE Y &EXIT 2 

&COMMAND XEDIT &FN FORTRAN A . 

&TYPE Do you wish to run the program again? (Y) 
&READ VAR &RESP4 

&IF &RESP4 EQ Y &GOTO -H 

&EXIT 2 

pees LOUR OUTPUT IS IN THE FILE &FN OUTPUT A 
&TYPE Do you wish to BROWSE your output? (Y) 
&READ VAR &RESP 

&IF &RESP EQ Y &COMMAND BROWSE &FN OUTPUT A 
Seer ee Print your output file? (Y) 

&READ VAR &RESP7 

2 ame EQ Y &COMMAND PRINT &FN OUTPUT A 
-REDO 

&TYPE Do you wish to KEDIT the program file? (Y/N) 
&READ VAR &RESPS5 

PiEeoRESPS EQ Y XEDIT &FN FORTRAN A 

&TYPE Do you wish to run the program again? (Y) 
&READ VAR &RESP6 

&RESP56 = &CONCAT OF &RESP5 &RESP6 

&IF &RESP56 EQ YY &GOTO -H 

&IF S&RESP6 EO Y &GOTO -RUN 

&EAIT 
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50 See need to compile your program ? (Y) 
§s FORTRAN COMPILER ENTERED. 12:11:11 
AAMAIN*X END OF COMPILATION 1 ****xx% 
AADEMAND** END OF COMPILATION 2 ***xxx 
ARARRIVE** END OF COMPILATION 3 ***xxx% 
AXASTART** END CF COMPILATION 4 ***xxx 


AROPTION** END OF COMPILATION 5 ***xxx 
VS FORTRANS COMET EERe 2. cue o. IZ2tiecs 


EXECUTION BEGINS... 


* Sue OPTION Bok DATA FILE (1 - 4) 
1 DATA1 2 DATAZ 3 = DATA3 4 = DATA4 


* ENTER RATE OF STOCKOUT COST TO UNIT PRICE ! ; 
1.25 0.91 = STOCKOUT COST PER UNIT IS 1252NGheUli@eeeRiGE 
a oaon eee STOCKOUT COST PER TIME IS 91% OF UNIT PRICE 


* ENTER OPTION FOR ORDER POLICY (1 = 2) ! 
1 = INDIVIDUAL ORDER POLICY 
2 = JOINT ORDER SOULE, 


CONSTANT LEADTIME 


NTER THE TYPE OF LEADTIME(1 - 2) ! 
= RANDOM (UNIFORM WITH MEAN) 


* ENTER SIMULATION EIMET YEARS ! 
Ql = SIMULATION LIMIT IS 1 YEAR 
OZ = SIMULATION Limi? £5 2 YeARs 


Ep 99 = SIMULATION LIMIT IS 99 YEARS 


* ENTER OPTION FOR START LEVEL eee 
i REGULAR(IND.ORDER POLICY: E.0.0 / JOINT ORDER POLICY: OD.UP QTY.) 
Z RANDOM BETWEEN (MUST ORDER LEVEL+1) AND (ORDER UP LEVEL) 


* ENTER OPTION FOR PRINTOUT(1 - 4) ! 
PRINT ALL ITEM 

PRINT ONE EVENT eo reM 

PRINT ONLY SUMMARY 

PRINT ONLY FINAL RESULT 


4 
YOUR OUTPUT 13> IN OTHE PILE ws Wi Oo re Ui 
Do you wish to BROWSE your output? (Y) 


bP WNrF 
Te | 


Print your output file? (Y) 

n 

Do you wish to XEDIT the program file? (Y/N) 
n 

Do you wish to run the program again? (Y) 


n 
R; T=4.63/5.82 12:11:54 
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** INPUT DATA : 


* USED DATA FILE NO. = 1 
ORPER FOLIGCYSOPITION = 1 
SIMIMULATION LIMIT(YRS) = 50 
START LEVEL OPTION = 2 
PRINTOUT OPTION = 4 
Pearse Tee ae=— 1 


a eS 


* ORDER POLICY OPTIONS 


Tee Da TO UAEVORDER POLICY 
Zs POSTE ORDER wrOoml CY 


* START LEVEL OPTIONS 


1 : REGULAR (IND. POLICY - E.0.Q ;: JOINT POLICY - OD.UP QTY.) 
2 : RANDOM (REORDER PONT + 1. ORDER-UP POINT ) 


* PRINTOUT OPTIONS 
Rewer LEM 
Zee ee ONLY EVENT ITEM 
See e eke ONLY SUMMARY 
4 : PRINT ONLY FINAL RESULT 
oeAb TIME TYeaS 


IGLOS Ue A) P29 BBN 8D 
2 : RANDOM LEADTIME (UNIFORM W/ MEAN) 


* SERVICE LEVEL: P2= 99 % SSA Ge boaAbe TIME: 0.083 YRS.( 1 MON. >} 
ITEM YR DEMAND Be MUST BUY CAN BUY ORD. UP START Exe eel . 
1 290.000 LS ZS DE ee) oy) 0308s 
Z 41.000 143 3 35 95 45 0.083 
3 77.000 109 7 48 93 38 8) 088 
- 22000 178 10 yet 150 60 O20 52 


* INPUT DATA FOR COSTS: 


Mbt UNIT PRICE H.RATE H.COST ORDER GOST B/O COST 
GROUP IND. Tie Tine 
HL 6.90 Qa20 es: 50.00 10.00) Rife O07 
Z tee 0) 6 218 0.24 50.00 10.00 0.48 0.01 
3 Smid O20 O26 50.00 (8) 101.0) ele 0.04 
“ 2.30 O20 0.46 50.00 10.00 OZ O02 
**X FINAL RESULTS (AVG.) : DURTIHG 50>" YEARS ARK 
ITEM # # DEMANDS # ORDERS AVG. UNIT YEAR AVG. B/O SERV ies 
ON HAND B/ O25 QTY. ee LEVEL (%) 
1 292.24 1.84 ood 0.025 2592 0.009 33.00 
2 41.58 G220 69 22s O08 0.34 0.006 eke le 
3 7 io O.72 54.064 0.004 0.46 OBE 10)S ee) 40 
x 120550 0.63 88.760 0.014 0.84 0.006 eho 5 2)(0. 
Met eaAvG.;- S3i.32 3.54 
AVG. NO. OF JOINT REP; 0.000 
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18243) 


PD LIN) 


* COMPUTED AVG. ANNUAL COSTS ($) 


HOLDING 


IOS e226 
LonoL2 
42.170 
40.830 


BACKORDER ORDER 
TY. slaligic 
5 OG. 036 110.400 
Oss Oe eZ 18.000 
O.7us8 0.003 43.200 
O62 S' 0.007 40.800 


TOTAL = 212.400 
STAND. DEVIATION = 
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TOTAL 


hes 
SA 
090 
409 


ose 
942 


kk INPUT DATA : 


USED DATA FILE NO. = 1 
ORDER PCLICY OPTION = 2 
SIMIMULATION LIMIT(YRS) = 
Sten Ie bevel, OPTION = 2 
PRINTOUT OPTION = 4 

Baap TIME Ivrea — 1 


A OO Ob Oo 


* ORDER POLICY OPTIONS 


i Lb eure ORDER POLICY 
2 OL IeRDER POLICY 


* START LEVEL OPTIONS 


fence ibs POLICY = E.0,0 = JOINT POLICY = OD. UP QTY.) 
2 : RANDOM (REORDER POINT + 1; ORDER-UP POINT) 


* PRINTOUT OPTIONS 
Peek aol LTE 
Zeno be ee VvENT . ITEM 
3 : PRINT ONLY SUMMARY 
4 : PRINT ONLY FINAL RESULT 
PEEAD TIME TYPES 


ie eCOostAN] LEADTIME 
2 : RANDOM LEADTIME (UNIFORM W/ MEAN) 


* SERVICE LEVEL: P2= 99 % AVG. HEAD Tite. 0.083 YRS.(¢ 1) MON.) 


ITEM YR DEMAND EQ MUST BUY CAN BUY ORD. UP START Ge 6 edhe 


l 290.000 15 25 ee ion eae 0.083 
2 41.000 143 -1 35 95 28 0.083 
3 G7. 00@ 109 5 48 93 32 0.083 
4 122.000 178 7 oa 150 50 0.083 
* INPUT DATA FOR COSTS: 
ITEM UNIT PRICE H.RATE H.COST ORDER COST B/O COST 
GROUP IND oer TIME 
iE 6.90 0.20 1.38 50.00 16.00 .76 Oro 
2 Pe 20 O20 O24 50.00 10.00 Onde Oa 
3 3.90 o720 0.758 50.00 10.00 1.56 0.04 
4 2.30 0.20 0.46 50.00 10.200 0-92 0) 4102 
*X FINAL RESULTS (AVG.) : DURING 50 YEARS ARK 
ITEM # # DEMANDS # ORDERS AVG. UNIT YEAR AVG. B/O SERVICE 
ON HAND B/O.S ge TIME LEVEL (%) 
iL 292.24 1.94 75.843 0.021 .56 0.007 99.12 
2 41.58 0.60 56.035 0.000 0.00 0.000 100.00 
3 7700 1.18 51.541 0.001 0,04) 0.001 99.95 
4 120.50 1.06 80.055 0.003 0,50" 0-007 99.75 
TOT AVG.: 531.32 4.78 
AVG. NO. OF JOINT REP; 1.620 
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ITEM 


iP WIN 


* COMFUTED AVG. ANNUAL COSTS (S$) 


HOLDING 


104.6€3 
13.448 
40.202 
SGro 2. 


ACKORDER ORDER 

BS Ge TIME 
066 02030 113.400 
0.000 0.000 6.000 
0.062 0.001 12-600 
O-276 0.001 PS c00 


TOTAL = 145.800 
STAND. DEVIATION = 
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TOTAL 


348. 
Zee 


i: 


i. 


i. 


13. 
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